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The  following  report  documents  the  design  study  of  the 
Air  Force  Institute  of  Technology  Graduate  Systems 
Engineering  Class  of  1985.  The  report  is  in  three  volumes. 
The  Executive  Summary  (Volume  I)  is  a  cursory  review  of  the 
study  and  is  meant  to  be  self-contained.  The  Final  Report 
(Volume  II)  and  the  Appendices  (Volume  III)  are  more 
detailed  and  should  be  read  together  for  completeness.  This 
study  explains  a  two-phase  methodology  we  developed  to 
permit  selection  of  an  optimal  military  satellite  servicing 
system.  The  work  was  conducted  from  December  1984  to 
December  1985.  The  original  project  concept  and  follow-on 
technical  support  was  provided  by  the  Rocket  Propulsion 
Laboratory  at  Edwards  AFB ,  California.  Additional  technical 
support  and  funding  was  provided  by  the  Office  for  Manned 
Spaceflight  (SD/YM)  and  the  Office  of  Plans  (SD/XR)  at  IJSAF 
Space  Division,  Los  Angeles  Air  Force  Station.  California. 
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Satisfaction  of  MP 1  for  Hierarchy  of  Objectives 


LEVEL  3 


Attributes:  Z1 

Z2 
Z3 
Z4 


Initial  Cost 
(Operating  Cost)/Time 
Rel iabi 1 ity  — 

(Mass  of  Payload  De 1 i vered ) /Time 


««******«  Hierarchy  Tree  of  Objectives  »»**»*** 


Satisfaction  of  mutual  preferential  independence  (MPI) 
among  the  attributes  of  an  objective  tree  is  a  necessary  and 
sufficient  condition  for  using  an  additive  value  function 
(Feldman  and  Rowell,  1985).  When  three  or  more  attributes 
are  involved,  satisfying  pairwise  preferential  independence 
(PPI)  is  equivalent  to  satisfying  MPI  (Keeney  and  Raiffa, 
1976:114).  PPI  and  MPI  are  defined  in  Section  3. 1.2. 2. 
Because  the  objective  tree  above  has  four  attributes,  and 


A-l 


since  satisfying  PPI  requires  fewer  steps,  satisfaction  of 
PPI  (which,  in  turn  satisfies  MPI)  is  demonstrated  below. 

In  order  to  satisfy  PPi,  the  analyst  or  decision  maker 
must  test  ail  possible  combinations  of  pairs  of  attributes 
for  at  least  two  levels  of  their  complementary  sets. 

This  requires 

n  1 

2  *  -  ( A .  1  ) 

2! (n-2) ! 

tests  (Helstrom,  1984:35).  In  this  case.  n  =  4,  and  a  total 
of  twelve  tests  are  required. 

Each  attribute  in  the  complementary  set  of  attributes 
should  be  set  at  a  given  level  while  the  preferences  between 
the  attribute  pairs  are  examined.  The  level  of  the  comple¬ 
mentary  set  should  then  be  changed  and  the  preference 
between  the  same  pair  of  attributes  should  then  be  re¬ 
examined.  If  the  preference  did  not  change  for  that  attri¬ 
bute  pair  when  its  complementary  set  changed  values,  that 
pair  of  attributes  can  be  considered  PPI  of  its  complemen¬ 
tary  set. 

It  is  recommended  that  the  two  levels  used  for  the  com¬ 
plementary  sets  be  the  least  acceptable  and  the  most  accept¬ 
able  values  for  those  attributes.  This  test  is,  of  course, 
a  subjective  one.  One  must  be  convinced  that  he  can  satisfy 
PPI  in  this  case  to  justify  using  a  linear  value  function. 
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The  tests  for  the  attributes  in  the  above  tree  are 


below.  where 

cost)/tirae,  "R" 

1 oad  de 1 i vered 

ATTRIBUTE 

PAIR 

”IC”  is  initial  cost, 

is  reliability,  and  ”MPD” 

to  orbit) /time . 

COMPLEMENTARY 

SET  AND  LEVEL 

”0C”  is  (operating 

is  (mass  of  pay- 

PREFERENCES 

CHANGES? 

UC.OC) 

{R,MPD> : 

lowest  values 

NO 

<IC,OC> 

<R,MPD>: 

highest  values 

NO 

{ I C . R> 

{OC . MPD> 

:  lowest  values 

NO 

{ I C , R) 

{OC , MPD) 

:  highest  values 

NO 

{  IC , MPD) 

<OC,R>: 

lowest  values 

NO 

{IC.MPD} 

<OC,R>  : 

highest  values 

NO 

<0C , R> 

{ IC ,MPD> 

lowest  values 

NO 

{OC,R> 

{IC.MPD) 

highest  values 

NO 

{OC , MPD> 

{ 1C , R>  : 

lowest  values 

NO 

<0C , MPD) 

{  IC  ,  R)  : 

highest  values 

NO 

{ R , MPD) 

{ IC , 0C>  : 

lowest  values 

NO 

<R.MPD> 

<IC.OC> : 

highest  values 

NO 

Therefore , 

PP1  and 

MPI  are  satisfied 

for  the  attributes 

in  this  hierarchy  tree,  and  a  linear  value  function  is  jus¬ 
tified. 


( 
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Appendix  B 


Detailed  Heirarchy  of  Objectives 
for  Selection  of  a 
Satellite  Servicing  System  (SSS) 


After  discussions  with  decision  makers  at  the  USAF 
Space  Division  (Green,  1985;  Lemon,  1985;  Sundberg,  1985; 
Wimberly.  1985;  Wittress,  1985;  Zerson,  1985)  a  comprehen¬ 
sive  hierarchy  of  objectives  was  developed,  consisting  of 
six  levels  and  40  objectives  and  subobjectives.  Figure  B.l 
shows  this  hierarchy  with  labeled  objectives  and  attributes. 


The  first  digit  of  the  numeric  label 
the  objective  within  the  hierarchy, 
used  to  identify  each  objective  with 
butes  are  found  in  the  circles,  and 


indicates  the 

level  of 

and 

the 

second 

digit  is 

i  n 

that 

1  eve  1 . 

Attr i- 

are 

identified 

by  a 

number  for  each  attribute,  preceeded  with  the  letter  ”A”  . 
The  following  is  a  description  of  each  objective  and  attri¬ 
bute  within  the  hierarchy. 


Descriptions  of  Objectives  Within  Figure  B.l 


1- 1  Military  Satellite  Servicing  System  Selection:  To 

select  the  satellite  servicing  system  which  best  meets 
the  objectives  presented  within  the  hierarchy  struc¬ 
ture  . 

2- 1  Satisfies  Congressional  Concerns:  The  system  is 

desired  to  meet  all  requirements  necessary  that  would 
impact  Congress’s  criteria  for  selecting  a  system. 

2-2  Mission  Accomplishment:  The  system  is  desired  to 
accomplish  its  mission  by  meeting  various  specifica¬ 
tions.  These  specifications  indicate  how  well  a  system 
could  accomplish  its  mission. 


2- 3  Utilization  of  Limited  Resources:  This  objective  is 

to  minimize  the  use  of  limited  resources  by  a  system. 
The  limited  resources  are  such  things  as  energy,  man¬ 
power,  money,  materials,  and  transportation  infrastruc¬ 
ture  . 

3- 1  Internationally  Politically  Stabilizing:  The  system 

is  not  desired  to  internationally  destabilize  relations 
that  exist  between  governments.  Such  things  as  a 
country’s  public  safety,  threat  of  attack,  or  breaking 
of  treaties  by  use  of  a  system  could  cause  political 
destabilization  of  international  relations. 

3~2  Technological  Risk:  This  objective  is  to  minimize 
risks  associated  with  selecting  systems  requiring  new, 
advanced  technology  developments  not  already  available. 

3-3  National  Prestige:  The  system  considered  is  desired 
to  favorably  impact  the  view  held  of  the  this  country 
by  others. 

3-4  Cost  Versus  Performance:  The  system  is  desired  to 
meet  its  performance  requirements  for  a  reasonable 
cost . 

3-5  Environmental  Impact:  The  system  should  not  adversely 
impact  the  environment. 

3-6  Deployable  Within  Constraints:  The  system  should  be 
able  to  be  implemented  within  apriori  time  require¬ 
ments  . 

3-7  Can  Supply  Satellite  Needs:  It  is  desired  that  the 
system  considered  be  able  to  meet  satellite  supply 
requ i rements  . 

3-8  Flexibility:  The  system  is  desired  to  be  flexible,  in 
that  it  is  able  to  operate  under  various  conditions. 

The  system  should  handle  anomolies  during  a  mission, 
and  be  able  to  respond  to  anomolies  by  priority. 

3-9  Survivability:  How  well  can  the  system  operate  under 
man-induced,  hostile  environment.  This  objective  is  to 
operate  "well”  under  hostile  environment. 

3-10  Availability:  is  the  system  available  based  on  its 
reliability  and  maintainability  characteristics?  The 
system  is  desired  to  be  available  over  as  much  time  as 
is  possible. 
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3-11  Transportation  Infrastructure:  The  system  being  con¬ 
sidered  is  desired  to  be  able  to  use  existing  transpor¬ 
tation  structures.  A  system  requiring  new  roads, 
canals,  ships,  airplanes,  to  transport  any  components 
would  not  be  as  preferred  as  a  system  not  requiring 
these . 

3-12  Energy:  This  objective  seeks  to  minimize  energy  use. 

3-13  Money:  This  objective  seeks  to  minimize  system  cost. 

3-14  Scarce  Strategic  Materials:  The  system  under  con¬ 
sideration  is  desired  to  use  a  minimum  amount  of  scarce 
strategic  materials  such  as  titanium,  nuclear  elements, 
or  chrome. 

4~1  Degree  of  Technology:  To  what  degree  does  the  system 
use  new  technology.  Does  the  system  utilize  advanced 
concepts  and  materials. 

4~2  Perception  of  System  Reliability:  The  system  is 

desired  to  be  perceived  as  being  very  reliable  in  its 
performance  of  activities. 

4-3  Americans  in  Space:  The  system  is  preferred  to  use 
increased  numbers  of  Americans  in  space  to  perform  its 
activities . 

4-4  National  Resource  Consum'pt ion :  The  system  under  con¬ 
sideration  is  desired  to  use  minimal  amounts  of  the 
countries  national  resources. 

4-5  Pollution:  This  objective  is  to  minimize  pollution 
created  by  a  system.  Such  things  as  water,  air,  and 
land  pollution  are  desired  to  be  minimized  in  a  system 
se 1 ect i on  . 

4-6  Public  Safety:  The  system  being  considered  is  desired 
to  not  harm  the  public  in  its  operation. 

4-7  Handle  Anomolies  on  Mission:  How  well  can  the  system 
under  consideration  respond  to  unplanned  events  while 
performing  a  mission?  The  more  a  system  is  able  to 
respond  to  such  unplanned  events  would  increase  its 
possible  selection. 

4~8  Respond  to  Anomolies  By  Priority:  How  well  can  the 
system  repond  to  unplanned  events  when  that  system  is 
not  performing  a  mission.  Greater  response  capability 
would  indicate  a  more  preferred  system. 
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4- 9  Reliable:  The  system  is  desired  to  be  reliable  in  its 

operat i on . 

4~10  Maintainable:  The  system  is  desired  to  be  maintain¬ 
able  to  keep  it  operating. 

4-11  System  Life  Cycle  Costs:  This  objective  is  to  minim¬ 
ize  the  life  cycle  costs  incurred  by  the  system  being 
considered . 

5~1  Spare  Availability:  Does  the  system  have  spare  parts 
available  to  repair  damaged  equipment  when  the  system 
is  on-roission?  A  more  preferred  system  would  have  such 
spare  parts  available. 

5- 2  Robotic  Level:  The  system  is  desired  to  have  sophis¬ 

ticated  robotic  technology  incorporated  in  its  design. 

5-3  fl  of  People  On-board  OSV:  The  system  being  considered 
is  desired  to  have  people  on  board  to  handle  unplanned 
events  occuring  on-mission. 

5~4  Spares  Availability:  Does  the  system  have  spare  parts 
available  to  repair  damaged  equipment  when  the  system 
is  on-mission:?  A  more  preferred  system  would  have 
such  spare  parts  available. 

5-b  Response  Time".  The  sytem  is  desired  to  be  able  to 
respond  quickly  to  anoraolTes  by  priority. 

5-7  Operations  Cost  Per  Time:  This  objective  is  to  minim¬ 
ize  operations  costs  per  time  for  a  system. 

5- 8  Amortized  Initial  Cost/Time".  This  objective  is  to 

maximize  the  system  life,  while  minimizing  the  initial 
cost . 

6- 1  On-orbit  Storage  Capacity:  The  system  is  desired  to 

have  storage  capacity  to  carry  available  spare  parts. 

6-2  R&D  Cost  Per  System  Life:  This  objective  is  to  mi  ,im- 
ize  the  research  and  development  costs  in  developing  a 
system  with  a  long  life.  The  system  is  desired  to  have 
minimal  research  and  development  costs,  and  a  "long” 


6-3  Production  Cost/Hardware  Life:  The  system  being  con¬ 
sidered  is  desired  to  have  small  production  costs  and  a 
"long”  hardware  life. 


Descriptions  of  the  Attributes  Within  Figure  B.l 


A1  Agrees  With  Policies  and  Treaties:  Does  the  system 
agree  with  international  policies  and  treaties  in  its 
operation? 

A2  Nan  On-board:  Is  the  system  manned  or  unmanned  when 

operating  in  space? 

A3  R&D  invested:  How  much  research  and  development  is 
required  for  the  system  to  be  developed? 

A4  Net  Savings  Due:  The  measure  of  savings  from  servicing 
the  satellite  with  the  system,  instead  ot  replacing  the 
satel 1 ite . 

A5  Time  From  IOC  to  FOC:  The  difference  in  time  from  Ini¬ 
tial  Operational  Capability  (IOC)  of  the  system,  till 
Full  Operational  Capability  (FOC)  of  the  system. 

A6  %  Satellite  Needs  Net  Per  Time:  The  percent  of  the 

required  satellite  needs  that  can  be  met  by  the  system. 

A7  Scenario:  Description  of  an  envrionment  that  the  sys¬ 
tem  is  operating  in  to  aid  in  measuring  survivability. 
The  system  could  be  operating  in  a  scenario  of  a  peace¬ 
time  environment,  a  nuclear  war  environment,  or  a  con¬ 
ventional  war  environment” 

A8  R&D  Invested:  How  much  research  and  development  is 
required  for  the  system  to  be  developed? 

A9  Fuel  Specific  Impulse:  The  specific  impulse  of  the 
fuel  for  various  components  of  the  system  requiring 
fuel  . 

A10  Structural  Nass  Ratio:  The  structural  mass  ratio  for 
the  various  components  of  the  system. 

All  OSV  Spare  Parts  Storage:  The  OSV  spare  parts  storage 
capacity  for  a  system  that  uses  OSVs. 

A12  Robotic  R&D  Costs:  The  OSV  costs  incurred  for  research 
and  development  in  automating  the  operation  of  the  OSV 
for  a  system. 

A13  #  Robots  On-board:  The  number  of  robots  required  on¬ 
board  an  OSV  in  the  performance  of  its  mission. 


r 
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A14  Scenario:  a  description  of  the  environment  within 

which  the  system  would  be  operating.  Some  examples  are 
a  peacetime  environment,  nuclear  war  environment,  or  a 
conventional  war  environment. 

A15  Depot  Storage  Capacity:  The  capacity  within  a  system 
to  store  necessary  spares,  in  space,  at  a  central  loca¬ 
tion. 

Alb  OSV  Payload  Capacity:  The  OSV  payload  capacity  which 
could  be  used  to  carry  necessary  spare  parts. 
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Computer  Program  for  Calculation  of  Figures  of  Merit 
Using  Linear  Value  Functions 


PROGRAM  FOM 

C  PROGRAM  TO  CALCULATE  FIGURE  OF  MERIT  FROM  NDSS 
C 

C  VARIABLE  DECLARATION 
INTEGER  NS.NP 

REAL  P( 100, 100) ,PMAX( 10) .PM IN ( 10) , V( 100, 100) ,QP ,TC, IC,0C.R,MPD 
REAL  INF0M1 , 1NF0M2.TF0M( 100) 

C  IDENTIFY  FILES  FOR  INPUT  NDSS  AND  OUTPUT  DATA 
OPEN  (UNIT=10,FlLE=’ndss) 

OPEN  (UNIT=  3,FILE= ’data’ ) 

OPEN  (UNIT=  4.FILE= ’ rank ’ ) 

OPEN  (UNIT=  7 , F I LE= ’ va  1  ’ ) 

REWIND  10 
REWIND  7 

READ  IN  «  OF  NONDOMINATED  SOLUTIONS  AND  U  OF  PERFORMANCE  INDICIES 
NS=69 
NP=4 


READ  IN  THE  PERFORMANCE  INDEX  MEASURES  FOR  EACH  SYSTEM  REPRESENTATION 
DO  50  J=1  ,NS 

READ  (10.1000)  (P(J. I ) . 1=1 ,NP)  _ 

1000  FORMAT  (4E17.9) 

50  CONTINUE 

WRITE  THE  PERFORMANCE  INDEX  MEASURES  READ  IN.  AND  THE  SYSTEM  « 

WRITE  (3.1100) 

1100  FORMAT  (’SYSTEM  «’ .9X, ’OPERATING  ’.9X.’ INITIAL  ’,8X. 

•RELIABILITY’ ,7X.  MASS  OF  ’) 

WRITE  (3.1101) 

1101  FORMAT  (20X. ’COST’ .13X, ’COST’ ,23X, ’OF  PAYLOAD) 

WRITE  (3.1102) 

1 102  FORMAT  ( 69X . ’ DEL I VERED ’ ) 

WRITE  (3.1105) 

105  FORMAT  ( 17X. ’DOLLARS/HR ’ ,9X,  DOLLARS’ . 1 IX. ’ PERCENT ’, 1 3X , ’KG’ ) 

WRITE  (3.*)  ” 

DO  60  J=1  .NS 

WRITE  (3.1110)  J , (P( J , I ) , 1=1 ,NP) 

110  FORMAT  (2X.I3.7X.4E17.9) 

60  CONTINUE 

FIND  THE  MAX  AND  MIN  VALUE  OF  EACH  PERFORMANCE  INDEX 
DO  70  I = 1 , NP 
DO  80  J=1 .NS 

PMAX(I)  =  AMAX1(PMAX( I ) ,P(J. I )) 

PMIN(l)  =  AMIN1 (PMIN( I ) ,P( J , I ) ) 

80  CONTINUE 
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70  CONTINUE 


C 

C  DETERMINE  LINEAR  VALUE  FUNCTION  VALUE  FOR  EACH  PERFORMANCE  INDEX 
C  IN  EVERY  SYSTEM 
DO  100  J=1 ,NS 
DO  90  1=1.2 

C  THIS  FUNCTION  GIVES  VALUES  THAT  MINIMIZE  COSTS 

V(J.I)  =  (PMAX(l)  -  P(J,I))  /  (PMAX(l)  -  PMlN(l)) 

90  CONTINUE 
C 

C  THIS  FUNCTION  GIVES  VALUES  THAT  MAXIMIZE  RELIABILITY  AND  MASS 
C  OF  PAYLOAD  DELIVERED 
DO  95  1=3.4 

V(J.I)  =  (P(J.I)  -  PMIN(D)  /  (PMAXU)  -  PMIN(D) 

95  CONTINUE 
100  CONTINUE 
C 

C  WRITE  VALUE  FUNCTION  VALUES  FOR  EACH  PI  IN  EVERY  SYSTEM 
WRITE  (3,*)  ’ ’ 

WRITE  (3.*)  ’LINEAR  VALUE  FUNCTION  VALUES' 

WRITE  (3.*)  ’ ’ 

WRITE  (3.1100) 

WRITE  (3.1101) 

WRITE  (3.1102) 

WRITE  (3.*)  ” 

DO  110  J=1.NS 

WRITE  (3.1210)  J, (V( J , I ) , 1=1 ,NP) 

WRITE  (7.1211)  (V(J.D.I-l.NP).J 

1210  FORMAT  (2X. 13.7X.4E17.9)  — 

1211  FORMAT  (4E13.7.3X, 13) 

110  CONTINUE 

C 

C  READ  IN  WEIGHTINGS  FOR  HIERARCHY  OF  OBJECTIVES  FROM  DECISION  MAKER 
C  PREFERENCES 

WRITE  (*.*)  ’WHAT  IS  THE  HIERARCHY  WEIGHTING  FOR’. 

* ’ OVERALL  PERFORMANCE?’ 

READ  («,*)  OP 

WRITE  (*,*)  ’OVERALL  PERFORMANCE  WEIGHTING  I S’. OP 
WRITE  (4.*)  'OVERALL  PERFORMANCE  WEIGHTING  I S’. OP 
WRITE  (*,*)  ’WHAT  IS  THE  HIERARCHY  WEIGHTING  FOR  TOTAL  COST?’ 
READ  (*.*)  TC 

WRITE  (*.*)  ’TOTAL  COST  WEIGHTING  IS’.TC 
WRITE  (4.*)  ’TOTAL  COST  WEIGHTING  IS' ,TC 

WRITE  <*.*)  'WHAT  IS  THE  HIERARCHY  WEIGHTING  FOR  INITIAL  COST?’ 
READ  (»,*)  1C 

WRITE  (*.*)  ’INITIAL  COST  WEIGHTING  IS’.IC 
WRITE  (4.*)  ’INITIAL  COST  WEIGHTING  IS’.IC 
WRITE  (*,*)  ’WHAT  IS  THE  HIERARCHY  WEIGHTING  FOR’. 

* ’OPERATING  COST/HR?’ 

READ  («.*)  OC 

WRITE  (*.*)  ’OPERATING  COST/HR  WEIGHTING  IS’.OC 
WRITE  (4.*)  OPERATING  COST/HR  WEIGHTING  IS’.OC 
WRITE  (*.*)  'WHAT  IS  THE  HIERARCHY  WEIGHTING  FOR  RELIABILITY?’ 
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READ  (.*.*)  R 

WRITE  (*.*)  ’RELIABILITY  WEIGHTING  IS’ .R 
WRITE  (4.*)  RELIABILITY  WEIGHTING  IS' .R 
WRITE  (*,*)  'WHAT  IS  THE  HIERARCHY  WEIGHTING  FOR'. 

* 'MASS  OF  PAYLOAD  DELIVERED?’ 

READ  (*.*)  MPD 

WRITE  (*.*)  MASS  OF  PAYLOAD  DELIVERED  WEIGHTING  IS’. MPD 
WRITE  (4.*)  MASS  OF  PAYLOAD  DELIVERED  WEIGHTING  IS’ .MPD 
C 

C  CALCULATE  THE  FIGURE  OF  MERIT  FOR  EACH  SYSTEM  REPRESENTATION 
C  USING  THE  WEIGHTINGS  AND  VALUES  APPLIED  TO  THE  HIERARCHY  TREE 
DO  130  1=1. NS 
INF0M1=0.0 
1NF0M2=0.0 

INF0M1  =  U0C*V(1 .1))-KIC*V(1 .2)))*TC 
1NF0M2=1 (R*V(  1 .3)  )-*-(MPD*V(  1,4))  )*QP 
TFOM( I )=INF0M1^INF0M2 

WRITE  THE  FIGURE  OF  MERIT  FOR  EACH  SYSTEM 
WRITE  (3,*)  ’ ’ 

WRITE  (3.*)  TFOM( I ) . ’ IS  THE  FIGURE  OF  MERIT  FOR  SYSTEM’. I 
WRITE  (4.*)  TF0M( 1 ) . ’ IS  THE  FIGURE  OF  MERIT  FOR  SYSTEM’.l 
130  CONTINUE 
END 
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Model  Equations 

This  appendix  presents  the  equations  used  to  model  the 
LG+OSV,  FHG+LG+OSV,  and  SB+FHG+LG+OSV  Satellite  Servicing 
Systems.  The  Performance  Index  equations  are  in  section  D.l 
and  the  State  and  Constraint  equations  common  to  the  three 
models  are  in  section  0.2.  Additional  equations  added  to 
this  common  set  for  the  LG+OSV,  FHG  +  LG+OSV ,  and 
SB+FHG+LG+OSV  models  are  in  sections  D.4,  D.5,  and  D.6 
respectively.  A  complete  descriptive  listing  of  all  vari¬ 
ables  (state,  intermediate,  and  exogenous)  is  in  Appendix  E 
and  the  intermediate  equations  are  in  Appendix  F. 


D.l  Perfomance  Index  Equations 

Z1  is  Operations  Costs 

Equation  D . 1 

Z1  =  1 04+ I  45+ 1 05+  1 46 

104  =  Launch  Site  Operations  Cost 

=  40*(U79/8760)*( 103**0.34)/ [ ( 1 / X335 ) **0 . 55 ] 
* ( X320 ) / X335 

105  =  LG  Fuel  Cost 
=  X326*U31 *X301 

146  =  0SV  Manned  Cost 
=  X500*X  555*1178 
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145  =  OSV  Fuel  Cost 
=  X52b*U22*X50 1 

Z2  is  Initial  Costs 

Equation  D.2 

Z2  =  107+ 106+ 1 47+ 1 48 

106  =  LG  R&D  Cost 

=  6500*1179*  (X325**0 .21 >*U37*U38*[EXP(2*X345)  1 

107  -  LG  Production  Cost 

=  X300* 12*(X325**0.56 ) *{ X300**  t  LN(U32)/LN(2.0) ]} 
*U79*[EXP(2*X345> ] 

147  =  OSV  R&D  Cost 

=  6500*U79*( 144**0.21 )*U29*U30*[ (U21 )**X555] 
*IEXP(2*X545) ] 

148  =  OSV  Production  Cost 

=  X500*(16.5)*< 144**0.56) 

*<X500**ILN(U23)/LN(2.0>  ]  >*1)79*  l  EXP  ( 2*X545 )  ] 


Z3  is  Rel iabi 1 ity 
Equation  D . 3 . 1 

Reliability  for  the  LG+0SV  and  FHG+LG+OSV  models. 
Z3  =  {1-[(1~X345 ) **X300 ]}*{l-[( 1 -X545) **X500 1 > 
Equation  D.3.2 

Reliability  for  the  SB+FHG+LG+OSV  model. 

Z3  =  <l-[ ( 1-X545)**X5001) 


Z4  is  Mass  Delivered  to  the  Satellites  per  Time  (Hr) 


Equation  D.4 
24  =  X50l*X510 

=  (#OSV  missions/time)*(OSV  payload) 

D.2  State  and  Constraint  Equations  Common  to  all  Models 


Equation  0.5 

Inequality  equation  placing  upperbounds  on  the  time 
between  LG  launches  due  to  maximum  time  people  can  remain  in 
space . 

1/X301  <-  U77 

(time  between  LG  launches) 

<-  (max  time  people  can  be  in  space) 


Equation  D.6  and  0.7 

Two  inequality  equations  placing  bounds  on  the  time 
between  LG  launches  from  a  specific  launch  site. 

X335  ->  1)48 
U49  ->  X335 


Equation  D.8 

Inequality  equation  placing  lowerbounds  on  the 
required  time  between  LG  launches  for  refurbishment  of  the 
LG  vehicle.  This  equation  assumes  the  launch  vehicle  must 
sit  on  the  ground  for  a  U40  time  period. 

X330  ->  U40 
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Equation  D.9 

Equality  equation  relating  the  LG  launch  rate  to  the 
launch  capability  of  the  launch  sites. 

X301  =  X320/X335 
#LG  launches  = 

(ttlaunch  s i tes ) / ( t i me  between  launches  form  a  site) 


Equation  D.10,  D.ll 

Equations  placing  upper  and  lower  bounds  on  the  LG 
vehicle's  structural  mass  ratio. 

LG  structural  mass  ratio  [128  =  X325/ (X326+X325) ]  becomes 

128  ->  IJ34 
128  <-  U36 


Equation  D . 12 

Inequality  equation  requiring  the  LG  to  have  one  or 
more  stages. 

X370  ->  1 


Equation  0.13,  D.14 

Inequality  equations  placing  upper  and  lower  bounds  on 
the  LG  reliability. 

X345  <-  1 
X345  >-  0 

(  - 
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Egua  t i on  U  .  L 5 


Equality  equation  equating  the  delta  velocity  the  LG 
must  experience  to  attain  orbit  to  the  thrust  the  mlg  must 
provide.  Equations  for  LG  delta  velocity  are  from  (Hill  and 
Peterson  ,  1970  )  . 

U33«U93*X370«LN{129/[ 1  28*  ( 129-1)+1]>  =  121 
(Vmlg)  =  (lsp  LG  f ue 1 ) * ( gra v i tat iona  1  constant) 


*(«  LG  stages )*LNl (LG  gross  lift  off  weight) 
/(empty  LG  mass+LG  payload)] 

129  =  [  (X325^X32b+X3l0)/X310]**(  1  /X370) 

121  =  SORT  [  1190  /  (  X360+U91 )]  +5486+ 1 152+1097 +0 
128  =  X325/(X325+X326) 


Equation  0.16 

Inequality  equation  placing 
mission  length  due  to  the  maximum 
space . 


upperbounds  on  the  0SV 
time  people  can  remain 


in 


X500/X501  <-  U77 


Equation  D .  17 

Inequality  equation  requiring  the  0SV  to  have  one  or 
more  waiting  orbits. 

X565  ->  1 


Equation  D.18,  0.19 

Equations  placing  upper  and  lower  bounds  on  the  CSV 
vehicle’s  structural  mass  ratio. 
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OSV  structural  mass  ratio  1142 


X525/ ( X526+X525 ) ] 


becomes 

142  ->  1)  12 
142  <-  U 1 3 


Equation  0.20,  D.21 

Inequality  equations  placing  upper  and  lower  bounds  on 
the  OSV  reliability. 

X545  <-  1 
X545  >-  0 


Equation  D . 22 

Equality  equation  relating  the  OSV  delta  velocity  for 
one  mission  to  the  OSV  payload  and  required  fuel  for  the 
mission. 

135  =  1 33»(  X560- 1 )  +•  1 32 

=  U24*U93*LN{[ I44  +  X526  +  X5101/E  1 44  +  X526  1 } 

(Vosv  required)  =  (Vosv  is  capable  of  providing) 

134  =  (Ul+U9l)*[(U2*X565-l)/( U2*X565 )]**(2./3.) 

133  =  2*ABS{SQRT  C  1)90/  (U1  +  U91 )  ] 

-SQRT[2*U90/(U1+U91)-U90/I34]> 

136  =  SQRT{2*ABS[  U90  /  (  U  l  +1)91  )  —  U 90 /  (  X360+U I  +2*1)91)1} 

137  =  SQRT{2*A  BSC  1)90/  ( X360+U9  l)-U90/(  X360+U 1 +2«U91  )  ]  } 
132  =  2*ABS{SQRT  [U90/(Ul+l)91)]-I  36+  1 37 

-SQRTC  U90/  (X 360+ 1)91  )  ]> 
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Equation  0.23 
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Equality  equation  relating  the  OSV  payload  to  the 
number  of  satellites  serviced  per  mission  and  the  mass 
delivered  to  each  satellite. 

X510  =  X561 *X560 
OSV  payload  mass  = 

(avg  mass  de  I  i  vered  )  *  (.9  satellites  serviced/mission) 


Equation  0.24 

Equality  equation  relating  the  number  of  OSV  missions 
per  time  to  the  number  of  OSV  vehicles  and  time  from  the 
start  of  one  mission  until  the  start  of  the  next  mission. 
X50 1  =  X500/( 151+150) 

80SV  missions/time  = 

( JIOSV  )  /  ( time  from  one  launch  untill  the  next) 

150  =  2*U92*{[(2*U91+Ul +X360 )/2. 1**1 .5}/SQRT(U90) 

+(X560-1 )*<2*U92«[ (U1+U91 )*«1 .51 
»  (  U2*X5b5~  1  .  )  /  [  l)2*SQRT  (  U90 )  ]}  +  X5bO*152 

151  =  X510/U25 


Equation  D . 25 

Equality  equation  relating  the  total  payload  mass  car¬ 
ried  by  OSV/time  to  the  total  mass  delivered  to 
sate  I  1 i tes / time. 

(X501  )»(X510)  -  (U2WX561  )/U6 
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Equation  D . 2b 


Equality  equation  relating  the  number  of  satellites 
serviced  by  the  OSV  to  the  total  number  of  satellites 
requiring  servicing. 

X501 *X560  =  U2/U6 

(#OSV  Mission/tiroe)*(  (ISat  serviced/OSV  Mission)  = 

(total  #sate 1 l ites ) / ( serv ice  interval) 


Equation  D . 27 

Inequality  equation  placing  lowerbounds  on  the  OSV-LG 
rendezvous  altitude,  FHG  and  LG  target  altitude,  and  SB 
location  altitude. 

X360  ->  U58 

*»»«**•****«***********«•***•*»*******•**«***•***■**■**■**■*«**•***•««* 

D . 3  LG+OSV  Model 

For  the  LG+OSV  model  the  following  equations  must  be 
added  to  the  group  of  common  equations. 


Equation  D . 28 

Equality  equation  requiring  the  LG  to  carry  enough 
payload  into  space  to  satisfy  OSV  and  satellite  mass 
requirements . 

(X301 )*(X310)  =  [ 141+X501-X510] 


(#LG  missions/time)*(LG  payload)  =  (OSV  needs) 

+(sate!lite  needs) 
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Equation  D . 29 


Inequality  equation  requiring  the  OSV  fleet  to  be  able 
to  hold  the  payload  from  one  LG. 

X500*(X510+X526+U19+X555*<U76+100/U77)/X501 )  ->  X310 
(#OSV)*(OSV  payload  +  OSV  needs)  ->  LG  payload 


Equation  D . 30 

Equality  equation  relating  the  LG  mission  rate  to  the 
number  of  LGs  and  the  time  between  launches  of  a  single  LG. 
X301  =  X300/  (  1 12+X33Q ) 

(8LG  launches/time)  = 

( SLG ’ s  ) / ( tota 1  time  between  launches  of  a  LG) 

113  =  {SORT  [U90/  (X360+U91 )  1+5486  +  1 152+1097}/  (2*1)93) 

114  =  X310/U25 

115  =  3. 

112  =  I  13+ I  14+ I  15 

**»#»»*«****«*«*»***»»»**#»*«**»***»*******«********«***»*■* 

D . 4  FHG+LG+OSV  Model 

For  the  FHG-LG-OMV  model  the  following  equations  must 
be  added  to  the  group  of  common  equations.  LG+OSV  equations 
D . 29  and  D . 30  are  also  required  in  the  FHG+LG+OSV  model. 


The  following  operations  costs  roust  be  added  to  the 
common  operating  costs.  Z1  becomes: 

Equation  0.31 

Z1  =  104+1 05 +145+ I46+I 10b+I 1 15+1 1 16 
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1106  =  Cost  of  FHG  Ground  Based  Energy/Time 
=  X 100* ( 2* I 105)*U111 

1115  =  FHG  Launch  Site  Ops  Cost/Time 
«  X120*U112*(X135**-2) 

1116  =  FHG  Vehicle  Costs/Time 

=  X100*<(1 .057*0.0624534*1103*1108) 

+(46330*1 (0.4536*1 113/9 )**0. 77 J )+(Ul 14*1 113)} 


Equation  D.32 

The  following  initial  costs  must  be  added  to  the  com¬ 
mon  initial  costs.  Z2  becomes: 

Z2  =  I07+I06+I47+I48+I118+I119+I120 

1118  =  FHG  Vehicle  R&D  Costs 

=  <13.235*0.0624534*1103*11081 
+ 1 74 14460+22600*(JT4536* 1113/9]} 

1119  =  Launch  Site  R&D  Costs  (launcher  and  power 

source  or  plant) 

=  U1 17*X360*( 1 103**3) 

1120  =  Launch  Site  Purchase  or  Fabrication  Cost 

-  X 1 20*{ 1 U 1 18* ( I  103**3) 1  +  IU1 19* (2*1 105 >**0.51} 


Equation  D.33 

Inequality  equation  for  the  minimum  amount  of  mass  the 
LG  must  carry  into  orbit.  The  right  hand  side  of  the  ine¬ 
quality  represents  the  mass  that  must  be  launched  by  a  iow-G 
launch  system. 

X30 1 *X3 10  ->  (U14*X561*U2/U6)+I40+(X555*X500*U76) 

D-10 


( U19*X501  ) 


(LG  carried  mass/time)  ->  (low-G  satelite  needs/time) 

+  (OMV  req  people  mass/tirae) 
+(OMV  life  support  mass/tine) 
+  (OMV  parts  needs/time) 


Equation  D . 34 

Inequality  equation  placing  a  lower  bound  on  the  time 
between  launches  from  a  particular  FHG  launch  site. 

XI 35  ->  U 1 20 


Equation  D.35 

Inequality  equation  placing  lower  bound  on  the  payload 
mass  of  the  FHG.  The  equation  guards  against  the  FHG  launch¬ 
ing  ridiculously  small  payloadsT 
X110  ->  U 123 


Equation  D . 36 

Equality  equation  defining  the  relationship  between 
the  required  FHG  mission/ launch  rate  and  the  required  time 
between  launches  at  a  site  and  the  number  of  launch  sites. 
X100  =  X120/X135 

#FHG  1 aunches/ time  =  (Dlaunch  sites) 

/time  between  launches  from  a  single  site 


Equation  D . 37  (revision  of  equation  D.28) 

Equality  equation  defining  the  relationship  between 
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the  mass  launched  into  space  by  the  two  launch  systems  and 


the  mass  required  in  space  by  the  satellites  and  the  OSV . 
(X301  )*(X310)+(X100)*(X110)  =  (  141+X501*X510) 

(<*LG  M  i  ss  i  on  / 1  i  me  )  *  (  LG  Payload) 

+  (8FHG  M i ss ions*Usab I e  FHG  Pay  1  oad  /  M  i  ss  i  on  ) 

=  (OMV  needs / 1 ime )  +  ( mas s  delivered  to  sate  II i tes/ t ime) 


D . 5  SB+FHG+LG+OSV  Model 

For  the  SBFHG-LG-OMV  model  the  fol  lowing  equations 
must  be  added  to  the  group  of  common  equations.  FHG+LG+OSV 
equations  D.34,  E.35  and  D . 36  are  also  required  in  the 
SB+FHG+LG+OSV  model . 

Equation  D .  38  (revision  of  D.T7  and  D.31) 

The  following  operations  costs  must  be  added  to  the 
common  operating  costs.  Z1  becomes: 

Z1  =  104 +105 +145 +146+11 06+1115+1116+11 30+ 1131 

1106  *  Cost  of  FHG  Ground  Based  Energy/Time 
»  X100*(2*I105)*U11I 

1115  =  FHG  Launch  Site  Ops  Cost/Time 
=  X 1 20* U 1 12*(X135**~2) 

1116  ”  FHG  Vehicle  Costs/Time 

=  X100*<(1 .057*0.0624534* I 1 03* 1 108) 

+  ( 46330* [ <0.4536*1 1 1 3/ 9 ) **0 . 77  I )  +  < U 1 1 4* l 113)} 
1130  =  SB  Manned  Cost/Time 
=  X400*X415»U78 
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1131  =  SB  Replenishment  Costs/Time  (for  f  ue  l -’-parts ) 

=  FHG5*U 1 34  +  1124*11139 

Equation  D . 39  (revision  of  equations  0.2  and  0.32) 

The  following  initial  costs  must  be  added  to  the  com¬ 
mon  initial  costs.  Z2  becomes: 

Z2  =  148+107+106+147+1118+1119+1120+1134+1135+1136 

1118  =  FHG  Vehicle  R&D  Costs 

=  [(3.235*0.0624534*1103*1108) 
+(7414460+22600*0.4536*1113/9)1 

1119  =  Launch  Site  R&D 

-  U117*X360*( 1 103**0.5) 

1120  =  Launch  Site  Purchase  or  Fabrication  Costs 

=  Xl20*{ [ U1 18* ( 1 103**0.5) ]  +  [U119*(2* 1 105 )**0. 51} 

1134  =  SB  R&D  Costs  “ 

=  U136*(X4 25**0. 5)+U135*(X415**0. 75) 

1135  =  SB  Production  or  Fabrication  Costs 
=  X400*U137*(X425**0.5) 

1136  =  SB  Development  Costs 
=  X400*U 1 38*X425 

Equation  D.40  (revision  of  equation  D.29) 

Inequality  equation  requiring  a  single  SB  to  be  capa¬ 
ble  of  holding  an  LG  entire  payload  and  the  required  SB 
saf  ty  1  eve  1 . 


X420  ->  X310+[U130*( 1 4 1 +  X5 1 0*X50l  +  I  124+1 I25+I 126+1 127) J 


Equation  D.41  (revision  of  equation  D.33) 

Inequality  equation  for  the  minimum  amount  of  mass  the 
LG  must  carry  into  orbit.  The  right  hand  side  of  the  ine¬ 
quality  represents  the  mass  that  must  be  launched  by  a  low-G 
Launch  system. 

X30 1  *X3 1 0  ->  (Ul4*X561*U2/U6)+I40+( X555*X500*U76 ) 
+(U19*X501)+I125+I 126+1127 

(LG  carried  mass/time)  ->  (low-G  satellite  needs/time) 

+  (OSV  required  people  mass/time) 
+(OSV  life  support  mass/time) 

+  (OSV  parts  needs/time)+(SB  needs) 


Equation  D.42 

Equality  equation  for  the  SB  structure  mass  to  SB  mass 
storage  capacity.  This  is  a  way^  of  insuring  the  SB  structure 
is  large  enough  to  hold  the  mass  the  SB  is  expected  to 
store . 

X425/X420  =  U 1 33 


Equation  D.43  (revision  of  equations  D.37  and  D.28) 

(X301 )*(X310)+(X100)*(X110)  =  ( I41+X501*X510) 

+1124+1125+1126+1127 

#l.G  mi  ss  ions  / 1  i  me  )  *  ( LG  payload) 

+(#FHG  mi ss i ons*usab 1 e  FHG  pay  1 oad / m i s s i on  ) 

=  (OMV  needs)+(mass  del  to  Sa t / t ime ) + ( SB  needs) 
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Equation  D.44  (revision  of  equations  D.30) 

Equality  equation  relating  the  LG  mission  rate  to  the 
number  of  LGs  and  the  time  between  launches  of  a  single  LG. 
X301  =  X300/ ( I 12+X330) 

(#LG  launches/time) 

=  ( #LG * s  ) / ( tota 1  time  between  launches  of  a  LG) 

113  =  {SQBT  [  1)90/  ( X360+U91 )  ]  +5486+ 1152+1 097}  /  (  2*1193  ) 

114  =  X310/U35 

115  =  3. 

112  =  1 13+ 1 14+ 1 1 5 
X301  =  X300/ ( I 12+X330) 

(#LG  launches/time) 

=  ( #LG ’ s ) / ( tota 1  time  between  launches  of  a  LG) 


D-15 


hysical  Model  Schematic 
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Appendix  E 


Model  Variables 

Low-G  State  Variables 
X300  -  If  of  LG  systems 

X301  -  ff  LG  missions/time  (mission  rate/hr) 

X310  -  Payload  mass  per  launch  [kg/ launch! 

X320  -  ff  of  launch  sites 

X325  -  Empty  vehicle  structure  mass  [kgj 
X326  -  Mass  of  propulsion  fuel  [kg] 

X330  -  Average  time  between  missions/LG  [hrs] 

(downtime  for  all  reasons) 

X335  -  Average  time  between  launches  from  a  specific 
launch  site  [hrs] 

X345  -  Reliability  of  LG  [0<X345<1] 

X360  -  OSV-LG  rendezvous  ahtitude 
X370  -  ff  LG  stages 

OSV  State  Variables 
X500  -  «  of  OSV s 

X501  -  ff  OSV  missions/time  (OSV  mission  rate) 

X510  -  Payload  mass  per  OSV  mission  [kg/missiom] 

X525  -  Mass  of  OSV  empty  [kg] 

(without  guidance  and  life  support  equip) 

X526  -  Mass  of  propulsion  fuel  req/mi ss ion( kg/mi ssion ] 
X545  -  Reliability  of  an  OSV  [0<X545<1] 

X555  -  OSV  crew  size  (ff  of  people) 

X560  -  ff  satellites  serviced/OSV  mission  [number/mission] 
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X561  -  Average  mass  delivered  to  a  sate  l  1  i te/ serv ice 
( mass / serv i ce ) 

X565  -  K  of  waiting  orbits  (n  in  the  orbital  equations) 

Fixed  High-G  State  Variables 
X100  -  H  of  FHG  launches/ time 
X110  -  Usable  payload  mass/launch 

(mass  usable  by  OSVs,  SBs ,  sate  1 1 i tes , etc ) 

X120  -  U  of  FHG  launch  sites 

X135  -  Time  between  launches  at  a  specific  site  Ihrsl 

SB  State  Variables 

X400  -  8  of  Space-bases  (SB) 

X415  -  8  of  people  to  operate  the  SB 
X420  -  SB  mass  storage  capacity 
X425  -  SB  structure  mass  — 

Exogenous  Variables 

The  exogenous  variables  are  input  variables  used  to  in¬ 
put  constants  into  the  model.  When  possible  the  exogenous 
variables  are  taken  from  some  reliable  source  document  or  an 
expert  in  the  area  the  exogenous  is  involved  with.  The  vari¬ 
ables  that  have  as  a  source  "group  est”  were  either  not 
found  in  the  literature  or  were  not  in  a  form  usable  in  the 
model.  In  these  cases  the  study  group  came  up  with  an  esti¬ 
mate  that  seemed  resonable  from  an  engineering  point  of 
view.  As  more  work  in  the  area  of  satellite  servicing  takes 
place  in  the  future  more  accurate  values  for  the  exogenous 
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can  easily  be  input  and  the  model  rerun  to  find  the  effect 
of  the  changes  . 

U1  -  Average  satellite  altitude  [kmj 

200<  HI  <  23000  nominal  =  800  kilometers 

Source,  group  est 
U2  -  #  of  satellites 

12<  U2  <  500  nominal  =  14 4  satellites 

Source,  group  est 

U6  -  Satellite  service  interval  [hrsl 

(time  between  services  to  the  same  satellite) 

8760<  U6  <  43800  nominal  *  2628013  yrsl 

Source.  General  Dynamics , 1983 
U12  0SV  mass  structure  ratio  lower  bound 
U12  -  0.1 

Source.  Hill  and  Pettrerson  ,  1970  :  328 
U13  0SV  mass  structure  ratio  upper  bound 
U 1 2  =  1.0 

Source.  Hill  and  Petterson , 1970 : 328 
U14  -  Percent  of  satellite  required  mass  which 
requires  low-G  launch 

0.0  <  U14  <  1.0  nominal  =  0.25  or  25% 

Source,  group  est. 

U16  -  Coef  for  computation  of  0SV  life  support  equiplkgl 
2432  =  U 1 6 

Source.  Bolen,  1980:Ch3 
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U 1 7  -  Coef  for  computation  of  OSV  life  support  equipment. 


[ kg/person ] 

305  =  U17 

Source.  Bolen,  1980:Ch3 
U18  -  Mass  of  OSV  guidance  equipment.  Ikgl 

185  <  U 1 8  <  277.5  nominal  =  200 

Source,  group  est. 

U19  -  OSV  mass  needs /mi ss i on : 

assume  all  needs  require  low-G  launch 

[OSV  parts,  expendables,  no  fuel  or  people!  [kg! 

100  =  U 19 

Source,  group  est. 

U21  -  Coef  effect  of  man  on  the  OSV  for  R&D  of  the  OSV 
0  <  U21  <  1  nominal  =  0.85 

Source,  group  est. 

1)22  -  COST/UNIT  OF  OSV  FUEL  [dollars/kg] 

$0.25  <  U22  <  $0.40  nominal  =  $0 . 32/kg 

Source.  Koelle,  1982-'  413 
Source.  Rehder,  1979 

U23  -  Cumulative  ave  learning  curve  slope  for  OSV 
0.95  -  U23 

Source.  Fong,  1981:  Ch4 
U24  -  OSV  ISP.  [hr  I 
U24  =0.14 

Source.  NASA.  1984d  :  42 
U25  -  LG  to  OSV  mass  transfer  rate 
(LG-0SV  and  FHG-LG-OSV  model) 


OSV  and  SB  mass  transfer  rate  (kg/time) 

450  <  025  <  900  kg/hr  nominal  =  800kg/hr 

Source.  Watkins,  1985 

026  -  OSV  to  Satellite  mass  transfer  rate  (kg/time) 

225  <  026  <  450  kg/hr  nominal  =  400kg/hr 

Source.  Watkins,  1985 

029  -  OSV  technical  development  correction  factor 
0.5  <  029  <  1.25 
new  system  --1.25 
technology  exists--0.8  to  1.0 
remake  existing  system-0.5  to  0.8 

Source.  Koelle,  1982:  405  nominal  =  1.25 

030  -  OSV  R&D  team  experience  factor 
0.7  <  030  <  1.3 

1.3*=no  experience  0.7  previous  experience 
Source.  Koelle,  1982:  405  nominal  *  1.3 

031  -  Cost/unit  of  LG  fuel  Idol  lars/kg] 

$0.25  <  031  <  $0.40  $/kg  nominal  =  $0 . 32/kg 

Source.  Koelle,  1982:  413 
Source.  Rehder,  1979 

032  -  Cumulative  average  learning  curve  slope  for  LG 
0.95  =  032 

Source.  Fong,  1981:  Ch4 
033  -  LG  ISP  I  hr  .  J 
0.11  =  033 

Source.  NASA,  1984d  :  54 


U34  LG  mass  structure  ratio  lower  bound 
U34  =0.1 

Source.  Hill  and  Petterson,  1970:  328 
U35  LG  to  SB  mass  transfer  rate  (kg/hr) 

450  <  U35  <  900  nominal  =  900kg/hr 

Source.  Watkins,  1985 
U36  LG  mass  structure  ratio  upper  bound 
U36  =  1  .0 

Source.  Hill  and  Petterson,  1970:  328 
U37  -  LG  technical  development  correction  factor 
0.5  <  1)29  <  1.25 
new  system  —  1 . 25 
technology  exists-~0.8  TO  1.0 
remake  existing  sys-0.5  TO  0.8 

Source.  Koelle,  1982:  405  nominal  =  1.25 

U38  -  LG  R&D  team  experience  factor 
0.7  <  U30  <  1.3 

1.3=no  experience-  0.7  previous  experience 
Source.  Koelle,  1982:  405  nominal  =  1.3 

1140  -  Min  LG  vehicle  required  downtiree/mission 
( land  to  launch) 

168  <  1)40  <  720hrs  tav  type  168  shuttle  type  720 
Source,  group  est.  nominal  =  400hrs 

U48  -  LG  launch  site  minimum  time  between  launches 
24  hrs  =  1)48 
Source,  group  est 
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U49  -  LG  launch  site  maximum  time  between  launches 
87b0hr  s  =  1)49 
Source,  group  est 

U58  -  Min  alt  LG  must  be  able  to  obtain 
1 85Km  =  U58 
Source,  group  est. 

1)75  -  SB  life  support  mass  requirements  per  person  time 
[ kg/ person*hr ]  (ie  a  i  r  ,  water  ,  food  etc) 

0.1  <  U75  <  1  nominal  =  0.20 

Source.  Guy,  1983  note'.  SB  will  recirculate 

(J76  -  Life  support  mass  requirements  per  person  time 
[ kg/ person»hr ]  ( ie  a ir , water , food  etc) 

0.1  <  U76  <  1  nominal  =  0.64 

Source.  Guy,  1983  note:  0SV  will  not  recirculate 
U77  -  People  time  in  space  before  rotation  [hrs] 

(radiation  exposer , physica 1 , 1 abor  contract,  etc) 
720  <  U77  <  4320  nominal  =  2160 

Source.  NASA.  1984b:  400 
U78  -  Cost/man-time 

$1000/man-hr  =  1)78 
Source,  group  est 

U79  -  Cost/man-year  for  cost  equations 

U79  =  $  1 25000/roan-year  1983  dollars 
r 

Source.  Koelle,  1982:  402 
U90  -  Earth  Gravitational  Constant(rou) 

=  5 . 16587 16e 12[km**3/hr**2l 
Source.  Bate,  1970:  429 
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U91  -  Earth  Radius  6378.145  km 
Source.  Bate,  1970:  429 
U92  -  PI  -  3.141592654 

Source.  Bate,  1970:  429 

U93  -  Earth  Gravitional  Constant  at  Earth  Surface(ge) 
=127137.6 

Source.  Hill  and  Petterson,  1970:  322 
Ulll  -  Cost  per  unit  of  earth  based  energy 
(initial  launch  energy] 

Ulll  =  $0 .00115/watt-hr 

Source.  one-half  Dayton  Ohio  Electrical 
rate  of  $2.30/kw-hr] 

U112  -  Coef  for  FHG  launch  site  operations  cost: 

U 1 12  =  $100000 
Source,  group  est. 

U114  -  Coef  for  cost/unit  of  FHG  apogee  fuel  burned 

$0.25  <  0114  <  $0.40  nominal  ==$0. 32/kg 

Source.  Koelle,  1982:  413 
Source.  Rehder,  1979 

U117  -  Coef  for  FHG  launch  site  R&D  costs - 
U 1 1 7  est.  *»  $100 
Source,  group  est. 

U118  -  Coef  for  FHG  launch  site  fabrication/purchase 
U 1 18  est  =  =  $i00 
Source,  group  est. 
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Cl  19 


U 120 


U  1 2 1 


U 123 


U  129 


U 1 30 


U  1 32 


Coef  for  FHG  power  >urce  fabrication/ purchase  est 
C  1  1 9  =  =  $100 
Source,  group  est. 

Min  time  between  launches  from  a  specific 
FHG  launch  site 
1)120  est.  ==  24hrs 
Source.  group  est. 

lsp  of  FHG  fuel  used  during  apogee  delta  vel 
[use  ave  chemical] 

U121  =0 . 09hr 
Source,  group  est.. 

Min  mass  of  payload  per  FHG  launch 
(usable  by  0MV,  satellites) 

U 123= 100kg 
Source,  group  est 

X  of  SB  structure  mass  for  fuel  calculation/time 
U 1 29  =  0.00000239 
Source.  NASA.  1984b:  465 

Coef  for  SB  safty  level  [length  of  time]  =lmonth 
U 1 30  =  720 

Source,  group  estimate 

X  of  SB  structure  mass  for  parts  caiculation/time 
this  represents  the  total  mass  of  the  SB 
replaced  every  lO.Oyrs 
U 1 32  =  0.000011415 


Source.  Space  station  configuration  discription 


U133  -  SB  structure  mass  to  mass  storage  capacity  ratio 
0.5  <  1)133  <  2  nominal  =  1.0 

Source.  NASA,  1984b:  45 

U134  -  Cost  per  unit  of  mass  for  SB  upkeep  (ie.  parts) 

U 1 34  =  $100. 00/kg 
Source,  group  estimate 

U135  -  Coef  for  SB  R&D  manned  influence 
U 130  =$500000.0 
Source,  group  est 

U13b  -  Coef  for  SB  R8.D  structure 
U 1 30  =$1000000.0 
Source,  group  est 

U137  -  Coef  for  SB  production  cost 
0130  =$1000000.0 
Source,  group  est 

1)138  -  Coef  for  SB  deployment  cost 
U 1 30  =  20000/kg 
Source,  group  est 

1)139  -  SB  station  keeping  fuel  cost/kg 

$0.25  <  1)139  <  $0.40  nominal  =  $0.32 


Source,  group  est 


Intermediate  Variables 

103  -  Launch  site  designed  capability 

(designed  launch  rate)  [ 1 aunches / t ime 1 

104  -  Launch  site  operational  cost/time  (for  each  site) 

105  -  LG  fuel  cost/time 
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LG  R&D  cost 


107  -  LG  production  cost 

112  -  LG  mission  time 

113  -  LG  -time  from  launch  to  rendevous  with  OSV  or  SB 

114  -  LG  docked  time 

(time  to  unload  LG  payload  to  OSV  or  SB) 

115  -  LG  -  return  to  earth  time 

121  -  LG  delta  velocity  per  mission 

128  -  Structural  mass  ratio  of  LG 

129  -  [LG  loaded  mass/payload  mas  s  ]  **  (  1  /  ft  LG  stages) 

132  -  OSV  delta  vel  for  resupply  (up  and  back)/0SV  mission 

133  -  OSV  delta  vel  in  waiting  orbit  (1  sat  to  1  sat) 

134  -  Semi-major  axis  of  elliptical  waiting  orbit 

135  -  OSV  total  delta  vel/OSV  mission 

136  -  Velocity  at  apogee  of  OSV  resupply  transfer  orbit 

137  -  Velocity  at  perigee  of  OSV  resupply  transfer  orbit 

140  -  OSV  req  people  mass  to  orbit/unit  time 

141  -  OSV  required  mass/time 

142  -  Structural  mass  ratio  of  OSV 

143  -  Mass  of  req  equipment  for  manned  life  support 

(nonrobotic  needs) 

144  -  Total  OSV  net  mass= 

structure  +  gu idance^- 1  i  f  e  support  equipment 

145  -  OSV  fuel  cost/time 

146  -  OSV  manned  cost/time 

147  -  OSV  R&D  cost 

148  -  OSV  production  cost 
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150  -  OSV  mission  length  <  time/OSV  mission 


=(time  to  travel  from  supply  orbit  to  satellites 
orbit  and  back)/0SV  mission 

+  l (total  time  to  travel  in  waiting  orbit  between 
satellites  serviced) 

+  (time  spent  servicing  satellites  on  location)! 

/ OSV  mission) 

151  -  OSV  resupply  time 

152  -  Time  to  service  one  satellite 

1100  =  FHG  vehicle  vel  leaving  the  ground  based  launcher 

1101  =  FHG  prelaunch  velocity  (earth’s  rotation  velocity) 
1103  =  Actual  mass  launched  from  FHG  earth  launcher 

1105  =  Min  ground  based  FHG  energy / 1 aunch 
to  provide  required  delta  velocity 
I  106  =  Cost  of  FHG  ground  based  energy/time 

1108  =  FHG  apogee  delta  velocity  required 

1109  =  Required  FHG  parking  orbit  velocity 

(so  FHG  will  stay  in  rendevous  orbit) 

1110  =  FHG  velocity  prior  to  appogee  burn 

1111  =  Parking  trendevousl  altitude  from  earth’s  center 
1113  =  Mass  of  fuel  required  at  the  FHG  apogee  burn 

1115  =  FHG  launch  site  operations  cost/time 

1116  =  FHG  vehicle  costs/time 

1118  =  FHG  vehicle  R&D  costs 

1119  =  FHG  launch  site  R&D 

(launcher  and  power  source  or  plant) 

1120  =  FHG  launch  site  purchase  or  fabrication  costs 
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1124  =  SB  fuel  required/time  (all  SB’s  in  operation) 

1125  =  SB  parts  mass  required/time 

(parts .nonfuel ,or  nonlife  support) 

1126  =  SB  life  support  mass  requ i red/ t i roe 

1127  =  SB  people  mass/time 

1130  =  SB  manned  cost/time 

1131  =  SB  rep  1 ent i shment  costs/time  ( f ue 1 +par ts+e tc  only) 

1134  =  SB  R&D  costs 

1135  =  SB  production-fabrication  cost 

1136  =  SB  deployment  costs 
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Appendix  F 


Intermediate  Variable  Equations 

The  intermediate  variables  used  in  the  model  are  for 
algebraic  simplification  only.  When  possible  a  name  is 
given  to  the  variable  as  an  aid  in  understanding  the  equa¬ 
tion  the  variable  represents.  When  a  literature  source  was 
identified  for  an  equation  the  source  is  provided.  In  other 
cases  when  no  literature  could  be  found  the  group  derived  an 
equation;  therefore,  as  future  work  and  research  continues 
the  equations  should  be  updated. 

103  =  LG  Launch  Site  Designed  Capability 

103  =  1.5/X335 

Source.  Group  Derived  (assumed  design  rate  =  1.5  req  rate] 

104  =  LG  Launch  Site  Operational  Cost 

104  =  40*(U79/8760)«( 103**0.34)/ ( ( 1/X335)**0.55] 

*  <  X320 ) / X335 

*  (40«(man  year  cost/hrs  in  yr) 

*( launch  site  designed  capaci ty*»0 . 34 ) 

*(l/(launch  site  launch  rate**0.55)J 

*(#  launch  sites)/(time  between  launches  at  a  site) 

Source .( Koe 1 1 e ,  1982:417)  cost  in  1983  dollars 

105  =  LG  Fuel  Cost/Time 

105  =  X326*U3 1 *X30 1 

-  fuel/mission  «<cost/unit  of  fuel)  *(#LG  missions) 

Source.  Group  Derived 
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10b  =  LG  R&D  Cost 


106  =  6500 *U79* 1 X325**0 .21  l*U37*U38*[exp(2*X345)  ] 

=  6500*U79* ( empty  mass  we i ght**0 . 2 1 ) *tech  factor 
•experience  factor 

•[cost  multiplier  for  increased  reliability! 
Source.  ( Koe 1 1 e , 1 982 : 407 )  cost  in  1983  dollars 

107  =  LG  Production  Cost 

107  =  X300« 12*1X325**0. 561*1 X300** { ln(U32)/ln(2.0)>l 
*U79*[exp(2*X345) 1 

=  (»LGs)*(  (  16.5+7.5)/2>*((X325)**0.56) 

*  ( (  #LG*»  (  1  nli32/  1  n2 )  )  *(cost  per  man-year) 

•[cost  multiplier  for  increased  reliability] 
Source.  Koelle,  1982:  409  cost  in  1983  dollars 

Source.  Rand,  1971:  Ch  5 

112  =  LG  Mission  Time 

112  =  113+114+15 
Source.  Group  Derived 

113  =  LG  Time  from  Launch  to  Rendezvous  with  0SV 

113  =  {SQRT  [  U90/  (  X 360 +  11 9 1  )  J  +5486+ 1 152+1097 }/(  2*1)93 ) 
Source.  Dept  of  AF,  1965:  2-46 

114  =  LG  Docked  Time  (time  to  offload  the  LG 

payload  to  0SV  or  SB) 

114  =  X310/L25  (for  LG+0SV  and  FHG+LG+0SV  model) 

114  =  X310/U35  (for  SB+FHG+LG+0SV  model) 

Source.  Group  Derived 
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115  =  LG  Return  Time  to  Earth 


115  =  3hrs 

Source.  Group  Derived 

121  =  LG  Delta  Velocity  per  Mission 

121  =  SQRT1U90/ (X360+U91 ) 1+5486+1 152+1097+0 

Source.  Dept  of  AF ,  1965:  2-46 

128  =  LG  Structural  Mass  Coefficient 

128  =  X325/ (X325+X326) 

Source.  Hill  and  Peterson,  1970:  328 

129  = 

129  =  X326+X31 0 )/X310J**( 1/X370) 

=  ILG  loaded  mass/payload  mass  1**(1/#LG  stages) 

Source.  Hill  and  Peterson,  1970:  328 

132  =  OSV  Delta  Velocity  (orbit  transition  from  resupply 
point  to  satellite  orbit  back  to  resupply  orbit) 

132  -  2*ABS(  SORT  (U90/(lil+ll91))-I 36+  I  37 -SORT  ( 1)90/  (  X360+U9 1  )  )  ) 
Source.  Appendix  G,  EQ(G.ll) 

133  =  OSV  Delta  Velocity  in  waiting  orbit  (1  sat  to  1  sat) 

133  =  2*ABS(  ( SORT (1190/  (U1  +  U91  )  )-SQRT(2*U90/  (U1  +  U91  )-U90/  134)  ) 
Source.  Appendix  G,  EQ(G.24) 

134  =  Semi-major  Axis  of  Elliptical  Waiting  Orbit 
134  -  (Ul+U91)»((U2*X565-l)/(U2*X565))**(2./3.  ) 

Source,  Appendix  G,  EQ(G.19) 
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135  =  OSV  Total  Delta  Velocity/OSV  mission 
135  =  1  33*  (  X560- 1  >-<-132 
Source.  Appendix  G,  EU(G.25) 

13b  =  Velocity  at  Apogee  of  OSV  Resupply  Transfer  Orbit 
1 3b  =  SQRT  (2*ABS(U90/  (U1+U91  )  -1190/  (X3b0-*-lil-"-2*li91  )  )  ) 
Source.  Appendix  G,  EQCG.8) 

137  =  Velocity  at  Perigee  of  OSV  Resupply  Transfer  Orbit 
137  =  SQRT(2*ABS(U90/(X3b0+U91)-U90/(X3b0  +  U1^2*li91))) 
Source.  Appendix  G,  EQ(G.b) 

140  =  OSV  -  People  Mass  to  Orbit/Unit  Time 

140  =  X500*X555* 1 00/ U77 
Source.  Group  Derived 

141  »  OSV  Required  Mass/Time 

141  =  l  (  X52b+U  19  )*X501  J  (  X555*U76*X500 )  +  1  40 
=  [ ( f ue 1 +par ts ) *#mi ss i ons 1 

-•■(life  support)-^ (people  mass) 

Source.  Group  Derived 

142  =  OSV  Structural  Mass  Ratio  (0.1<142<1) 

1 42=*X525/  (  X52b  +  X525 ) 

Source.  Hill  and  Peterson,  1970:  328 

143  =  Mass  of  Life  Support  Equipment  for  Manned  OSV 
143  =  U 1 6  +  U 1 7*X555 

Source.  Bolen,  1980:  Ch  3 
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144  =  OSV  Single  Vehicle  Mass  (dry  weight) 

144  =  X525+  U 1 8+  143 

=  structure+gu i dance+ 1 i f e  support  equipment 
Source.  Group  Derived 

145  =  OSV  Fuel  Cost/Time 

145  =  X526*U22*X501 

=  fuel  burned / m i ss i on* ( cost / un i t  of  fuel) 

*(#OSV  missions) 

Source.  Group  Derived 

14b  =  OSV  Manned  Cost/Time 

146  =  X500*X555*li78 

=  ( #OSV ' s )*(OSV  crew  s i ze ) * ( cost / mant i me ) 

Source.  Group  Derived 

147  =  OSV  R&D  Cost 

147  =  6500*1)79* l  144**0.21  ]*U29*(J30*[  (1)21  )**X555] 

*[exp( 2*X545 )  1 
=  6500*cost  per  man-year 

*(total  empty  mass  we i ght**0 . 2 1 ) 

*  tech  factor*  experience  factor 

*  [multiplier  for  addition  of  man  in  the  OSV 1 

*  [cost  multiplier  for  increased  reliability! 

Source.  ( Koe 1 le , 1982 : 407 )  in  1983  dollars 

148  =  OSV  Production  Cost 

148  =  X500*[ 16. 51*[ 144**0. 5bJ*[X500*»( ln(U23)/ ln(2.0) ) 1 
*U79*lexp(2*X545) ! 
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=  (»OSV)*(( 16. 5)*( 144**0.56) 

*(cumlative  ave  learncurve)*  cost  per  man  year 
*[cost  multiplier  for  more  reliability! 

Source.  Koelle.  1982:  407  in  1983  dollars 

Source.  Rand,  1971:  Ch  5 

150  =  OSV  Mission  Length 

=  (time  to  travel  from  supply  orbit  to  satelite 
orbit  and  back)/0SV  mission 
+l(total  time  to  travel  in  waiting  orbit 
between  satellites  serviced) 

■•■(total  time  spent  servicing  satellites 
on  location)]  /OSV  mission} 

150  =  2*U92*(((2*l)91+U  1 +X360  )/2.  )  **  l  .5)  /StlRK  U90  )  ) 

+  (  X560-  1  )*(2*U92*(  ( (j  1  09 1  )««1 . 5 ) * (  U2*X565~  1  .  ) 

/ ( U2*SQRT ( U90 ) ) ) *X560*152 
Source.  Appendix  G,  ECKG.35) 

151  =  OSV  Resupply  Time 

151  =  X510/U25 
Source.  Group  Derived 

152  =  Time  to  Service  One  Satellite 
152  -  X561/U26 

=  mass  del i vered 

/OSV  to  Satellite  mass  transfer  rate 
Source.  Group  Derived 
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1100  =  Post  Launch  FHG  Vehicle  Velocity 

(ie.  velocity  leaving  the  ground  based  launcher) 

1100  =  SORT  {2*ABS[  090/ 1)9 1 -1)90/ (2*091+X3b0)  ]> 

(note:  this  represents  the  speed  the  mass  (vehicle) 
would  be  traveling  to  be  at  the  perogee  of  an 
el i pt ica 1  orbit . ) 

Source.  Appendix  G,  EQ(G.b) 

1101  =  FHG  launcher  velocity. 

(speed  launcher  is  traveling  due  the  earths 
rotation,  note  our  launcher  is  considered  to 
be  physically  located  on  the  equator.) 

1101  =  U9l*2*U92/24 
=  U91  *1)92/  12 

=  ( rad i us ) * ( angu 1 ar  vel) 

=  (earth  radius  at  equator)*2 
*pi /( rotation  period  earth) 

Source.  Appendix  G,  EQ(G.5) 

1103  =  Total  FHG  Vehicle  launch  mass 
=  XllG  f  10*1113/9 

=  (mass  payload)  +  (10/9)*(mass  fuel) 

(note:  mass  structure  =  mass  fuel/9| 

Source.  Group  Derived 
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1105  3  Ground  Based  FHG  Energy  for  Launch 

(note:  for  simplicity  we  will  consider  only  the 
change  in  kinetic  energy  of  the  mass  launched 
from  earth) 

3  4 .629*0.5*1 103*( 1 100-1101 )  **2 
3  2.3145*1  103*( I  100 -1 101  )**2 
=  [ one-ha  1 f*mass*de 1 ta  vel**21 

•conversion  factor  to  get  1105  in  watt-hrs  (4.629) 
Source.  Group  Derived 

1106  =  Cost  of  FHG  Ground  Based  Energy/Time 

(note:  consider  an  electrical  source) 

(computed  based  on  twice  the  minimum 
required  energy  from  1105) 

=  X 100*  ( 2*  1  105  )  *U1 1 1  =cost  / 1  iroe 
Source.  Group  Derived 

1108  =  FHG  Apogee  Delta  Velocity  Required 
=  1109-1110 

Source.  Appendix  G,  EQ(G.10) 

1109  =  Required  FHG  Parking  Orbit  Velocity 

(so  FHG  will  stay  in  rendevous  orbit) 

=  SQRKU90/  1111} 

Source.  Appendix  G,  EQ(G.9) 

1110  3  FHG  Velocity  Prior  to  Apogee  Burn 

(ie  velocity  from  initial  launch) 
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SQRT(2*ABS[ U90/ 1111 


U90/(U91  +  1111)]} 


Source.  Appendix  G,  EQ(G.8) 

1111  =  Altitude  of  Parking  Orbit 

(rendevous  altitude  from  the  center  of  the  earth) 

=  U91  +  X360 

Source.  Group  Derived 

1113  =  Mass  of  Fuel  Required  at  FHG  Apogee  Burn 
=  X110*<  [exp( 1108/ (U121*U93) ) ] 

/  [  10/9  -  exp(I108/((J121*U93))/9  ]} 

(note:  this  equation  takes  into  consideration  the 
mass  of  the  unloaded  vehicle  by  assuming  the  mass 
structure  ratio  will  be  0.1.  This  assumption  is 
implemented  in  the  10/9  and  the  exp  (  )/9. 

Source.  Group  Derived 

1115  -  FHG  Launch  Site  Operations  Cost/Time 
=  X 120  *  U 1 12  *  (X135**-2) 

Source.  Group  Derived 

1116  *  FH  Vehicle  Costs/Time 

(note:  assume  FHG  Vehicle  is  expendable) 

=  X100»((  1  . 057  *0.0624534*1  103*1108) 

+(46330*1 (0.4536*1 1 13/9)**0.77J )  +  (UU4M113)} 

=  #FHG  1 aunches / t i me 

*(  engine  cost+p 1 atf orm  cost  +  fuel  cost} 

=  8FHG  1 aunches / t i me  *{(cost  coef) 

*(coversion  kg-km/hr  to  Ibf-sec) 
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*  (mass  to  delta  vel)  *(delt  vel  reg)l 

+  (cost  coef ) * ( kg  to  lfm  conv  *platform  mass)**. 77] 
+  [FHG  fuel  cost!) 

Source,  bong,  1981'.  Ch  5 

1118  =  FHG  Vehicle  K&U  Costs 

=  <[3.235*0.0624534*1 103*1 108) 

+  17414460  f  22600*0.4536*1113/9]) 

=  (engine  platform  R&D) 

Source.  Fong,  1981:  Ch  5 

1119  =  Launch  Site  R&D 

=  U1 17*X3b0*( 1 103**0 .5) 

Source.  Group  Derived 

1120  =  Launch  Site  Purchase  or  Fabrication 
=  X120*([U118*( 1103**0.5)  ] 

*  [UU9*(2*I105)**0.5J> 

Source.  Group  Derived 

1124  =  SB  Fuel  Mass  Required  per  Time 

*  X400*X425*ll  129  (note:  could  be  high-G  launched) 
Source.  Group  Derived. 

1125  =  SB  Mass  Requ i rements /Time  (require  low-g  launch) 

«  X400*X425*1J  132 

*=  [*SB’s]*(SB  structure  mass] 

*[%  SB  structure  mass  needing  replacement  per  time! 
Source.  Group  Derived. 
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1126  =  SB  Life  Support  Mass  Requ i red / T i me 
=  X400*U75*X415 

=  [  8  SB  ’  s  1  *1  ttpeop  le/SBl*Cmass/ person) 

(note:  assume  SB’s  recirculate  life  support  mass 
and  all  life  mass  is  low-G  launched) 

Source.  Group  Derived. 

1127  =  SB  People  Mass/Time 
=  X400*X415*100/U77 

Source.  Group  Derived 

1130  =  SB  Manned  Cost/Time 
=  X40Q*X4 15*078 

Source.  Group  Derived 

1131  =  SB  Rep  1 ent i shroent  Costs/Time  (for  fuel  +  parts  ) 

=  1 125*1)  134  +  1  124»U139 

*  l rep  1 en i shroent  mass / 1 i me  1  * { cost  of  mass] 

Source.  Group  Derived 

1134  =  SB  R&D  Costs 

=  U 1 36*{ X425**0 . 5}  +  IJ 1  35* { X4 1 5**0 . 75 > 

Source.  Group  Derived 

1135  =  SB  Product ion-Fabr i ca t i on  Cost 
=  X400*U  137*1  X 425**0 . 5  ) 

Source.  Group  Derived 

1136  =*  SB  Deployment  Costs 

*  X 400*  1!  1 38* X 4 25 

*  1  MSB ' s ]*[ launch  cost/kg]*  (SB  structure  mass) 
Source.  Group  Derived 
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Appendix  t> 

Derivation  of  OSV  Orhital  Mechanics  Equations 

The  size  of  an  OSV  is  determined  primarily  bv  two  fac¬ 
tors:  the  cargo  mass  it  must  deliver  to  the  satellite(s) 

and  the  fuel  mass  the  OSV  must  deliver  and  consume  to  accom¬ 
plish  the  mission.  Here  we  present  tne  development  of  the 
equations  used  to  model  the  dynamics  of  the  OSV.  Since  the 
SSS  models  require  OSV  quantities  on  a  per  mission  basis, 
the  following  development  is  for  a  single  OiV  mission, 
first  we  will  begin  by  defining  an  OSV  mission  and  present 
the  assumptions  made  to  simplify  the  form  of  the  equations. 

Then  the  equations  will  be  developed  for  the  OSV  fuel  con¬ 
sumption.  elliptical  orbits,  transfer  orbit,  waiting  orbit, 
and  mission  length/time  of  flight  (TOE).  Finally,  a  brief 

summary  of  the  primary  equations  is  given,  included  is  a  i 

cross-reference  table  to  the  intermediate  variable  equations 
(Appendix  F)  and  associated  model  equations  (Appendix  U). 

li .  1  OSV  Mission  Definition 

For  development  of  the  Satellite  Servicing  System  iSSS) 
model,  the  mission  of  an  OSV  is  defined  as:  ( 1 )  loading  on 
supplies  while  in  some  circular  parking  orbit.  iZ)  transfer¬ 
ring  from  the  parking  orbit  to  the  servicing  (satellite) 
orbit.  (3)  delivering  an  average  amount  of  mass  to  each 
satellite,  and  (4)  returning  empty  to  the  parking  orbit. 
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G . 2  Simplifying  Assumptions 


> 


To  develop  the  analytical  equations  that  describe  the 
OSV  mission,  the  following  simplifications  and  assumptions 
were  made : 

1.  100  percent  of  an  OSV  payload  capacity  is  distributed 
per  mission. 

2.  The  parking  orbit  is  the  same  at  the  beginning  and  end 
of  a  mission,  however,  it  is  tree  to  be  chosen  within 
the  mode  l  . 

3.  The  satellite  and  parking  orbits  are  circular. 

4.  There  are  no  perturbations  of  orbits  (transfer,  park¬ 
ing.  or  servicing). 

5.  The  effect  on  satellite  orbits  from  oblateness  of  the 
earth  is  neglected. 

b.  The  drag  of  the  atmosphere  and  solar  wind  on  OSV  is 
neg 1 ected . 

7.  Station  keeping  maneuvers  by  OSV  are  not  considered. 

8.  Rendezvous  maneuvers  are  neglected;  the  rendezvous  is 
assumed  to  occur  i mmed iately  at  intercept  of  transfer 
orbit  and  target  orbit. 

9.  Satellites  are  equally  spaced  within  the  orbit  plane. 

10.  No  priority  servicing  of  satellites  within  constella¬ 
tion.  Therefore,  the  OSV  will  travel  from  satellite  to 
satellite  in  a  consecutive  order. 

11.  The  maneuver  used  to  travel  between  satellites  within 
the  same  constellation  will  be  the  same  between  any  two 
sate l 1 i tes  . 

12.  The  requirements  of  satellites  in  a  constellation  will 
be  averaged  to  allow  the  same  amount  of  mass  delivered 
to  each.  Or  looked  at  another  wav.  satellites  within  a 
constellation  are  identical  and  therefore  require  the 
same  amount  of  mass  when  serviced. 

13.  In  order  to  keep  the  mission  time  length  short,  the  UbV 
propulsion  system  is  assumed  to  be  chemical. 
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G . 3  OSV  Fuel  Consumption 


Since  the  type  of  propulsion  system  is  a  chemical 
(solid  or  liquid),  the  fuel  mass  used  can  be  calculated  from 
the  Rocket  equation  (Hill,  1970): 


AVqsv  =  ge* 1 sp* 1 n ( M i /Mo  ) 


(G.  1  ) 


where 


ge  =  gravitational  acceleration  at  surface 
of  the  earth  (9.807  m/sec^) 

Isp  =  specific  impulse  of  fuel(sec) 

Mi  =  initial  total  vehicle  mass(kg) 

Mo  =  final  total  vehicle  mass(kg) 


AVqSv  =  magnitude  of  total  change  in  velocity 

of  OSV  to  travel  between  two  points(m/sec) 


Since  the  fuel  mass  needed  for  the  model  is  that  per  OSV 
mission,  then  the  final  and  initial  masses  of  the  OSV  are 
def i ned  as  : 


Mo 

=  Ms+Mp 

(G. 

.2) 

Mi 

=  Ms+Mp+Mf 

Mi 

=  Mo+Mf 

(G. 

.3) 

where 


Ms  =  OSV  structura 1 ( dry  )  mass  (kg) 

Mp  =  total  pay load( suppl y )  mass  delivered 
to  satellites  per  OSV  mission  (kg) 

Mf  =  fuel  mass  (kg) 
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Substituting  Eq  (G.3)  into  Eq  (G.l)  and  rearranging,  the 
fuel  mass  used  per  (JSV  mission  in  terms  of  the  USV  final 
(empty)  total  mass  and  the  payload  mass  is- 

Mf  =  Mo{exp[AVosv/(ge*lsp)  ]-l>  (G.4) 

The  fuel  mass  per  mission  calculated  by  using  Eq  (G.4) 
is  conservatively  high.  The  equation  is  derived  assuming 
that  the  payload  mass  is  removed  from  the  USV  all  at  once  at 
the  end  of  the  mission,  rather  than  dropping  mass  off  at 
each  satellite.  If  this  staging  process  were  to  be 
modelled,  an  iterative  application  of  Eq  (G.4)  would  be  used 
with  new  values  for  Mi  and  Mo  determined  at  each  start/stop 
point  of  the  OSV  mission.  Then  the  fuel  mass  used  for  each 
leg  of  the  mission  is  totalled  to  give  the  total  fuel  mass. 

A  EUKTRAN  program  called  VORBCS.F  is  developed  to  calculate 
the  fuel  mass  taking  into  account  staging  (see  Appendix  li  J  . 
Table  G.l  shows  a  comparison  of  the  fuel  mass  calculations 
from  this  program  with  those  using  the  approximation  of  Eq 
(G.4).  The  values  in  the  first  four  columns  describe  the 
mission  and  are  used  to  calculate  AVqsv.  These  equations 
will  be  described  later.  However,  notice  that  the  estimated 
fuel  mass  is  indeed  conservatively  high  and  the  difference 
between  the  fuel  mass  values  does  increase  with  increasing 
average  mass  per  satellite  and  number  of  satellites  ser¬ 
viced.  This  is  expected  because  the  effect  of  staging  saves 
more  fuel  mass  for  a  given  £-vosv  as  the  number  of  off¬ 
loading  events  increases.  In  Chapter  V  of  this  study 
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Table 


Number  or  waiting  orbits  (n)  set  =  1 
Altitude  of  parking  orbit  set  =  185.2 
OSV  structural  mass  (Ms)  set  =  2000  kg 
Isp  set  =  450  sec 


is  addressed  the  effect  of  using  Eq  (.6.4;  as  an  approxima¬ 
tion  by  examining  the  sensitivity  of  the  final  results  to 
changes  in  totai  fuel  mass  consumed. 

In  order  to  use  Eq  (G.4>,  the  AV qsv  must  be  calculated. 
This  is  dependant  on  how  an  OSV  moves  in  space.  We  treated 
the  motion  of  the  OSV  relative  to  the  center  of  the  earth  as 
a  classical  two  body  problem  in  orbital  mechanics.  The  mis¬ 
sion  dynamics  of  the  OSV  were  broken  up  into  two  distinct 
phases  :( 1  )  transfer  to  and  from  the  service  orbit  from  and 
to  the  parking  orbit;  and  (2)  transfer  between  satellites  of 
the  same  orbit,  i.e.  i ntersate 1 1 i te  maneuvers.  But  before  we 
discuss  the  AV  calculations  for  the  above  two  phases  we  must 
define  some  elements  of  the  two-body  problem;  namely  the 
elliptical  orbit  elements. 

G . 4  Elements  of  Elliptical  Orbits. 

For  a  small  body  (relative  to  earth’s  mass)  in  free 
flight  orbit  near  the  earth,  the  two  body  problem  predicts 
that  the  trajectory  of  the  body  will  be  that  of  an  ellipse 
with  the  center  of  the  earth  at  one  focus  (Bate.  1971:30). 
Figure  G.l  shows  the  elliptical  orbital  elements. 
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Figure  G.l  Elliptical  Orbit  Elements  (Dept,  of  AF  .  iyb5:2~lo) 
where 

a  =  semi  major  axis 
b  ~  semi  minor  axis 

c  *  distance  from  origin  to  either  focus 
ra  *  radius  of  apogee 
=  radius  of  perigee 

Notice  that  circular  orbits  are  a  special  case  of  elliptical 

orbits  when  r  =r  =r. 

a  p 

Since  the  OSV  is  assumed  to  have  chemical  propulsion, 
the  forces  used  to  change  its  velocity  will  be  considered 
impulsive.  Therefore,  after  an  instantaneous  change  in  velo¬ 
city,  the  OSV  will  be  in  free  flight  about  the  earth  and 
consequently  in  an  elliptical  orbit. 

G . 5  Transfer  Orbit 

There  are  many  ways  an  OSV  could  change  from  one  circu¬ 
lar  orbit  to  another  via  an  elliptical  transfer  orbit.  The 
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choice  depends  on  a  trade~off  between  and  transfer  time. 
The  quicker  one  wishes  to  change  between  orbits  the  more  AV 
(and  therefore  fuel)  the  vehicle  must  expend.  As  a  first 
modeling  attempt  for  this  study,  minimizing  fuel  consumption 
was  felt  more  important  than  minimizing  transfer  time. 


Therefore,  for  the  case  of  transfering  between  two  co- 
planar  circlar  orbits  a  Hohmann  transfer  was  modeled.  This 
roanuever  is  shown  in  figure  G.2  and  requires  the  minimum 
total  AV  (Bate,  1971:163). 


Figure  G.2  Hohmann  Transfer  between  Coplanar  Circular  Orbits 


The  Hohmann  transfer  trajectory  is  an  ellipse  such  that 
the  perigee  is  tangent  to  the  inner  circular  orbit  (rp=rpQ) 
and  the  apogee  is  tangent  to  the  outer  circular  orbit 
(r  =r__) .  Two  impulses  (or  burns)  are  required  to  complete 

8  SO 

the  maneuver. 
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Assume  the  transfer  is  from  the  parking  (inner)  orbit 


to  the  servicing  (outer)  orbit.  Figure  G.2  depicts  the 
orbital  velocities  relative  to  the  center  of  the  earth.  In 
order  for  the  OSV  to  initiate  travel  along  the  transfer 
ellipse  at  perigee  it  must  increase  its  speed  by  AVp.  The 
OSV  velocity  in  the  circular  parking  orbit  is 

V  =  (u/r  )l/2  (G  .5 

po  t'  po' 

where 

o  2 

=  Gravitational  parameter,  3.986032  x  10  ki  sec 
fpQ  =  Radius  of  parking  orbit 

and  the  OSV  velocity  in  the  elliptical  transfer  orbit  at 
perigee  would  be: 

V  =  {2[i|/r  -n/  ( r  +r  )  ] }  ^^  (G.6 

trp  x  ,r  po  r  po  so/J/ 

where 

rgQ  =  Radius  of  the  service  orbit. 

Therefore,  for  the  first  impulse  to  get  the  OSV  into  the 
transfer  orbit: 


IVtp  -VBnl 
trp  po 


( G  .  7  ) 


The  second  impulse  is  used  to  recircularize  the  OSV 
orbit  into  the  service  orbit  when  the  OSV  is  at  apogee.  The 
velocity  of  the  OSV  at  apogee  is  given  by: 
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and  the  velocity  of  the  circular  service  orbit  is 


therefore 


V 

so 


1/2 


so 


(G.  10) 


AV. 


=  IVso-Vtra' 


So  the  total  change  in  velocity  of  the  OSV  to  transfer 

one  way  between  two  orbits  is  AV  +^V  •  Therefore,  to  go  up 

pa 

to  the  service  orbit,  recircularize,  go  back  to  the  parking 
orbit,  and  recircularize,  the  total  change  in  velocity  is: 


AV 


resupply  *  21AVa+AVpl 


21Vso-Vtra+Vtrp-Vpo' 


(G  .  1 1 ) 


The  above  equation  is  good  only  for  orbits  in  the  same 
plane,  i.e.  zero  inclination  difference  between  them.  For 
orbits  at  different  inclinations  the  choice  of  transfer 
orbit  is  not  as  simple. 


The  three  roost  common  transfers  between  inclined  circu¬ 
lar  orbits  are  the  bi-elliptic  transfer  (Figure  G.3),  the 
"modified”  Hohmann  and  the  Hohmann  transfer  with  plane 
change  (Figure  G.4). 
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Figure  G.3  Bi-elliptic  transfer  orbit  (Dept  of  AF  .  19bo) 

For  orbit  transfer  without  rendezvous  the  Hohmann 
transfer  with  plane  change  uses  the  minimum  AV .  but  since 
the  OSV  must  begin  and  end  at  the  line  of  nodes,  a  satellite 
must  be  at  the  precise  angle  so  that  both  vehicles  arrive  at 
the  other  node  simultaneously.  The  bi-elliptic  transfer  and 
’’modified”  Hohmann  don’t  have  this  problem  but  do  use  more 
AV .  The  phasing  time  increases  as  the  ratio  of  the  final  to 
initial  orbit  radius  approaches  unity  (Baker.  1965:12).  How¬ 
ever.  it  is  assumed  that  the  ratio  of  parking  orbit  radius 
to  service  orbit  radius  will  not  be  near  unity  tor 
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Figure  G  .  4  "Modified”  Hohnann  and  Hohmann  transfer 
with  plane  change.  (Dept  of  AF.  1965) 

our  SSS  model.  In  addition,  the  first  satellite  to  be  ser¬ 
viced  in  an  orbit  is  arbitrary,  so  the  phasing  time  is  again 
reduced.  Since  it  was  felt  that  minimizing  fuel  consumption 
(or  )  was  very  important,  we  chose  the  Hohmann  transfer 
with  plane  change  to  model  transfers  between  inclined  circu¬ 
lar  orbits.  If  OSV  time  of  flight  becomes  more  critical, 
then  a  choice  between  the  other  two  techniques  is  recom¬ 
mended.  See  (Baker.  I9t>5)  for  a  detailed  discussion  in  this 
area  . 
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The  Hohmann  transfer  with  plane  change  is  a  two  burn 
maneuver  optimally  splitting  the  inclination  change  at  each 
burn  to  achieve  minimum  AV  .  The  first  burn  usually  includes 
an  inclination  change  of  less  than  5  degrees  with  the  rest 
occurring  at  the  second  burn. 

To  calculate  the  required  AV  at  each  burn  refer  to  fig¬ 
ure  G.5.  This  shows  the  velocity  vectors  at  a  point  of 
i ntersect ionn  between  two  orbits,  and  tnf  the  AV  required  to 
change  from  orbit  1  to  orbit  2.  Once  the  angle  0  is  known, 
the  AV  is  calculated  using  the  law  of  cosines  as  shown. 


Figure  G.5  Velocity  Vector  Diagram  (Dept  of  AF.  1 985 '.  2-40 ) 


Each  burn  of  the  Hohmann  transfer  with  plane  change 
combines  a  plane  change  maneuver  and  an  altitude  change 
maneuver  into  one.  Figure  G.b  shows  the  vector  diagram  of 
such  a  maneuver  at  the  intersection  of  a  transfer  ellipse 
(begun  at  100  nm)  and  a  final  circular  orbit  at  1 500nm  with 
an  inclination  difference  of  10  degrees. 

At  the  line  of  nodes,  a  in  Figure  G.b  is  the  inclina¬ 
tion  difference  between  the  two  orbits.  Let  0  be  the  amount 
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Figure  G.6  Combined  Maneuver  Vector  Diagram 
(Dept  of  AF,  1985:2-41) 

of  plane  change  (less  than  5  degrees)  at  the  first  burn  and 
•be  the  total  plane  change.  Then  the  change  in  velocity  at 
the  first  burn  is 

AV1  *  '''2p„-V2trp-2V[,0Vtrt>  coS(a>]1/2  <G.12) 

and  at  the  second  burn 

AV2  =  [v2Vtra+v2so-2^raVso  cos («-/?)]  1 ' 2  (G.13 

where  VpQ ,  Vtpp  ,  Vtra,  and  Vso  are  defined  by  Eqs  (G.5). 
(G.6),  ( G . 8 ) ,  and  (G.9)  respectively.  Therefore,  the  total 
change  in  velocity  between  inclined  orbits  for  up,  recircu¬ 
larize,  back  and  recircularize  is: 


AVresupply  =  2IAV1+£v2(  (G.14) 

Notice  when  9  is  equal  to  zero  this  equation  reduces  to  the 
co-planar  case  Eq  (G.ll).  One  of  these  two  equations 
((G.ll)  or  (G.14))  is  used  to  model  the  AV  used  for  phase  1 
maneuvers.  Now  let  us  turn  our  attention  to  modeling  phase 
2  maneuvers:  intersatellite  transfers. 


G-14 


Once  the  OSV  is  in  the  circular  servicing  orbit  and 
finishes  servicing  the  satellite  it  rendezvoused  with,  it 
must  move  on  to  the  next  satellite.  This  is  accomplished 
through  an  epoch  change  maneuver  (Dept  of  AF  .  19b5).  This 
maneuver  uses  the  same  principle  as  the  coplanar  Hohmann 
transfer.  An  impulse  is  used  to  place  the  USV  into  an 
elliptical  orbit.  Since  the  period  of  the  ellipse  will  be 
different  than  the  period  of  the  circular  satellite  orbit, 
the  OSV  and  satellites  will  move  with  respect  to  one 
another.  Figure  6.7  shows  the  case  for  an  ellipse  outside 
the  circular  orbit. 


Figure  G.7  Epoch  Change  Using  Outside  Transfer  Orbit 

( Dept  of  AF .  I9b5 ) 


The  size  of  the  ellipse  is  chosen  so  a  satellite  and 
the  QSV  rendezvous  after  an  integer  number  of  elliptical 
orbits  by  the  OSV .  Theoretically,  very  little  impulse  can  be 
used  on  the  first  burn  and  it  would  take  an  infinite  number 
of  revolutions  of  the  OSV  before  it  could  dock  with  a  satel¬ 
lite.  Conversely,  a  large  initial  impulse  could  be  used  to 


achieve  rendezvous  after  one  elliptical  orbit.  Here  again, 
the  trade-off  between  AV  and  time  must  be  made.  However,  the 
minimum  AV  requires  a  stiff  penalty  in  infinite  time. 
Therefore  a  decision  to  model  the  trade-off  was  made. 

An  epoch  change  can  occur  using  an  ellipse  inside  or 
outside  the  service  orbit.  For  every  exterior  ellipse  there 
is  an  interior  ellipse  that  requires  the  same  totalAV.  The 
choise  of  which  ellipse  is  determined  by  whether  the  OSV 
must  "catch  up  to”  or  "wait  for”  the  next  satellite.  Since 
there  is  no  predetermined  order  in  our  mode  1  we  chose  to 
investigate  using  an  interior  ellipse  because  it  has  the 
smaller  orbital  period.  This  choice  requires  that  precau¬ 
tions  be  taken  to  avoid  choosing  an  ellipse  which  could 
result  in  the  OSV  impacting  tke  earth.  Figure  G.8  shows  the 
geometries  used  for  the  equation  derivations. 


Figure  G.8  Waiting  Orbit  Geometry 
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In  order  to  assure  rendezvous  after  an  integer  number  of  USV 


waiting  orbits  the  following  must  hold: 

At  =  n(TSO-TwoJ  (G<15 

where 

At  =  Cha racterst i c  time  between  satelites 
(see  Figure  G  .  7  ) 

n  =  Interger  number  of  elliptical  waiting  orbits 

T  =  Papi  rk-f  n n p  corv  i  rp  orbit 

SO  - 

Twq  =  Period  of  one  waiting  orbit 

For  equally  spaced  satellites  At=Tso/N,  where  N  is  the 
number  of  satellites  in  the  orbit.  Using  this  in  Eq  (G.14) 
and  rearranging,  we  get: 


Two  =  Tsor (Nn-l)/NnJ  (G. lb) 

Which  clearly  shows  the  dependance  of  the  size  of  the  wait¬ 
ing  orbit  to  n;  the  number  of  waiting  orbits  before  rendez¬ 
vous.  The  period  of  each  orbit  is  related  to  the  orbit 
geometry  as  follows: 

TSo  -  2WrJ,2l^  ce.l7> 

tko  •  frri3l2iS}T 

where 

a  =  Semi-major  axis  of  waiting  orbit 


Using  these  two  expressions  in  Eq  (G.ib)  and  solving 
for  a  we  get : 
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a  =  rgQ( c  N  n - 1  >/Nni2/3 

s  i  nee 


(G  .  19) 


2a 


r 

P 


rVr 


2a~r 


P 

so 


r  +r 
so  p 


iG .  20) 


This  can  be  used  to  prevent  letting  the  Ubv  crash  into  the 
earth  by  placing  bounds  on  the  adius  at  perigee: 


rearth  T  185  km  <  rp  <  rSQ  1G.21) 

To  calculate  the  AV  to  insert  the  OSV  into  an  ellipti¬ 
cal  waiting  orbit  and  recircularize  it  after  n  ellipses 
refer  to  Figure  G.9 


Figure  G.9  Waiting  Orbit  Velocities 
The  velocity  of  the  service  orbit  is  repeated  for  conveni¬ 
ence  : 
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so 


(Jl/rso> 


1/2 


(  G  .  9  ) 


V 


The  velocity  of  the  waiting  orbit  at  apogee  is 


^woa  ~  {2  l^i  /  rgo  ^ r  sot  rp  ^ 


i/2 


(G.22) 


Therefore,  the  change  in  velocity  required  to  place  the  OSV 
into  a  waiting  orbit  is: 


AV  =  v  -v 

a  so  woa 


( G . 23 ) 


Noting  that  rso+tp  =  2a  the  total  change  in  velocity  to 


get  into  a  waiting  orbit  and  get  out  of  it  is: 


1/2, 


(G.24) 


So  ^Vwo  represents  the  total  AV  to  travel  from  one  satellite 
to  the  next.  Therefore,  if  Y  represents  the  number  of 
satellites  serviced  in  the  orbiX.  and  Y<N,  then  the  total  AV 
per  mission  is: 

AVoav  '  AVreSupply*a-1,AV»o  '■*■**> 

where  ^vresuppjy  is  given  by  either  Eq  (G.ll)  or  Eq  (G.14) 
and  AVwq  is  given  by  Eq  (.G.24).  These  equations  represent 
an  OSV  mission  profile  of  using  a  Hohmann  transfer  with 
plane  change  between  orbits  and  an  epoch  change  using  an 
inner  ellipse  for  intraorbit  changes  to  service  Y  satel¬ 
lites.  Using  this  profile  an  average  mission  duration  can 
be  estimated  for  the  SSS  model. 
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G  •  7  Mission  Length/Time  of  Flight  (TUF) 


Te  time  length  per  OSV  mission  can  be  conveniently 
broken  up  into  three  parts:  (1)  phas i ng/ a  1 i gn i ng  time,  (2) 
transfer  time  to  and  from  service  orbit  and  (3)  time  spent 
travelling  between  and  servicing  ¥  satellites  in  the  con- 
ste 1 1  at i on  . 

The  phasing  time  is  the  time  after  receiving  supplies 
that  the  OSV  must  loiter  in  the  parking  orbit  until  it  can 
initiate  orbit  transfer.  The  initial  burn  must  occur  at  a 
time  such  that  rendezvous  is  possible  at  the  second  burn 
(see  Figure  6.4).  This  launch  window  for  a  single  intercep¬ 
tor  vehicle  (OSV)  attempting  rendezvous  with  a  single  target 
vehicle  (satellite)  occurs  once  every  synodic  period  (Bate. 
1971:367).  The  synodic  period_for  earth  orbits  is  defined 
as  : 

Ts  =  2  TT/  I  wj-w^l  ( G  .  2b  ) 

where 

Tg  =  Synodic  period 

=  Angular  velocity  of  orbit  1  relative  to  the 
center  of  the  Earth 

w2  =  Angular  velocity  of  orbit  2  relative  to  the 
center  of  the  Earth 

But  the  angular  velocity  of  an  orbit  is  27T  divided  by  the 
orbit  period.  Therefore,  the  synodic  period  can  be  written 
in  terms  of  the  parking  orbit  and  service  orbit  periods. 
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1/  I  1/T 


V  (a  •  2.1  ) 


po  '^so1 


However,  in  our  model  every  satellite  is  a  candidate  target 
vehicle,  so  the  time  between  launch  windows  is: 


t  .  =  T  /N 

maxphase  s 


=  1/ 1 N/T  -  M/T  1 
po  1  ' ‘so 


(.  G  .  28  > 


where 


N  =  Number  of  satellites  in  the  constellation 


This  time  represents  the  longest  wait  before  next 
transfer  initiation.  It  was  assumed  a  launch  opportunity 
had  just  been  missed  at  the  completion  of  loading  the  OSV . 
Since  the  OSV  loading  completion  time  is  random,  we  chose  to 
use  the  mean  time  between  launch  windows  tor  our  static 
model.  Therefore,  the  average  phasing  time  per  OSV  mission 


is  one  half  of  t 


maxphase 


'avgphase 


1/ !2N/Tpo  -  2N/Tsol 


(G.291 


Substituting  T=27T rJ/ for  each  orbit  period  and  rear¬ 
ranging  the  above  equation  becomes: 


'avgphase 


■  (2Jr/77r)((rpo.rso)3,2/2Nlrso3';!-rp<)3/2.)  (S.30) 


This  equation  will  be  used  to  determine  the  time  the  OSV 
loiters  in  the  parking  orbit.  Now  the  transfer  times  will 
be  calculated. 


Since  all  maneuvers  are  impulsive  the  transfer  time 
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from  the  parking  orbit  to  the  service  orbit  and  back  again 
is  just  the  period  of  the  transfer  ellipse: 


Notice  this 


by  (r  -rr 

so  1  po 


Ttr  =  C27T/7)D[(rso+rpo)/2J3/2  <G.31> 

is  Eg  (G.18)  with  the  semi-major  axis  a  replaced 
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The  time  the  OSV  spends  in  the  constellation  can  be 
broken  into  time  spent  servicing  Y  satellites  and  the  time 
to  travel  between  them.  If  we  assume  some  average  time  to 
service  a  satellite,  say  Q,  then  the  total  time  to  service  V 
sate  1 1 i tes  i s : 


t 


servicing 


YU 


( G . 32 ) 


The  time  to  travel  between  Y  satellites  is  the  number 
of  trips  (Y-l)  multiplied  by  the  time  per  trip;  which  is 
simply  the  period  of  the  waiting  orbit.  Substituting  Eq 
(G.19),  the  semi-major  axis  of  the  waiting  orbit,  into  Eq 
(G.18),  the  period  of  the  waiting  orbit,  and  multiplying  the 
result  by  (Y-l)  we  get: 


Sntersat  =  <  2  7T /V^)  r  SQ3/ 2 U  Nn- 1 )  /  N  J  ( Y- 1 )  (G.33) 

Therefore  the  time  to  travel  between  and  service  Y  satel- 
1 ites  i s  : 


■  (2  7T^)rso:i/2aNn-l./»J(V-l)  ,  YQ  IG.34) 
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Finally  the  time  of  flight  or  mission  length  of  an  GSV 
is  the  sum  of  Eqs  (G.30),  (G.31)  and  (G.34): 

TOF  -  (27r/^T)l<r1)0.rso)3/2/(2«llrso:J^-rpo^2)l 

*  (2  7r/J}T)t(rsoTrp<,)/213/2 

♦  (27T/>/u)r  3/2[(Nn'l)/Nl(V-l)  +  YQ  (G.35) 

s  o 

G . 8  Summary  of  Equations 

We  present  here  for  convenience  those  primary  equations 


developed  in  the  previous  sections  that  areused  in  the  SSS 
model  to  represent  the  orbital  mechanics  of  the  OSV .  See 
appropriate  development  paragraphs  for  explanation  of  terms. 

Fuel  Consumed  by  an  OSV  per  Mission: 


Mf  -  Mo<exp[AVosv/(ge*lsp) J-1>-YP 


(  G  .  4  ) 


where 

P  =  average  delivered  payload  per  satellite 
Delta  Velocity  Used  by  OSV  per  Mission'- 


AV 


OSV 


■  ^resupply4  ^vwo 


( (.  .  2b  > 


where  A V 

wo 

and  ( G  .  2  i  ) 
the  Earth, 
transfers  a 


is  defined  by  Eq  (G.24)  and 
are  used  to  prevent  the  OSV 

AVresupply  is  defined  by  Eq 
nd  Eq  (G.14)  for  non~coplana 


Eqs  tG . 19) 
from  crash 
(li.lU  to 
r  transfer 


.  (G.20). 

ing  onto 
r  coplanar 

s  . 
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Time  Length  ot  an  USV  Mission: 

*  '■2lr,/Fn(rso*’-voi/^3ri 

+  (27r/^I)rso3/2t(hin-i)/N)(Y-l>  •*■  YU  (6.35) 

The  following  table  identifies  those  equations  derived 
in  this  appendix  which  are  used  as  intermediate  variable 
equations.  Also  identified  are  the  numbers  of  the  model 
equations  that  use  the  intermediate  variables. 


Table  6.2 

Equation  Cross-Reference 


Appendix  G 
Orbital 
Mechanics 
Equation  8 

Appendix  F 

I ntermediate 
Variable 
Equation  8 

Appendix  D 
Associated 

Mode  1 

Equation  8 

G .  4 

— 

D.22* 

G .  5 

1  101 

D  .  3 1 .  32.  38.  39 

G.b 

137.  100 

D.22.  31.  32.  38.  39 

G  .  7 

136,  110 

D.31.  32.  38.  39 

G .  9 

l  109 

D . 3 1 ,  32.  38.  39 

G.  10 

1  108 

D.31.  32.  38.  39 

G.  11 

132 

D.22 

G.  19 

134 

D.22 

G.25 

135 

D.22 

G .  35 

150 

D.22 

*  The  right  hand  side  of  D.22  is  a  rearranged  version  of  6.4 
in  terms  of  vehicle  structural  mass  ratio  (.  Ms  /  ( Mf +  Ms  )  )  la¬ 
beled  ! 4?  This  ratio  is  bounded  to  range  within  technolog¬ 
ically  achievable  limits. 
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Appendix  H 


Computer  Program  VURBCS.F  and  Sample  Output 


****************  double  precision  version  ********************* 

***********************  vorbcs.f  ****************************** 
*************************************************************** 
*************************************************************** 
*********  THIS  PROGRAM  IS  AN  ORBITAL  TRANSFER  MODEL.  ********** 
****  original:  27  MAY  1985  ;  this  version:  18  JUNE  1985  **** 

*********  CAPT  R  LIEBER  GSE-85D  ********** 

*************************************************************** 
*************************************************************** 

*  THIS  MODEL  CALCULATES  THE  DELTA  V  AND  TIME  OF  FLIGHT  FOR  * 

*  AN  OSV  TRAVELING  BETWEEN  A  RESUPPLY  ORBIT.  AND  A  SERVICING  * 

*  ORBIT  DETERMINED  BY  THE  ALTITUDE  AND  INCLINATION  OF  THE  * 

*  SATELLITES.  IT  IS  INTERACTIVE  AND  ALLOWS  THE  USER  TO  SPEC-  * 

*  IFY  THE  ALTITUDE  AND  INCLINATION  OF  THE  SATELLITE  CONSTEL-  * 

*  LATION.  THE  ALTITUDE  AND  INCLINATION  OF  THE  RESUPPLY  ORBIT.* 

*  THE  TOTAL  8  OF  SATELLITES  IN  THE  SERVICE  ORBIT.  THE  s  OF  * 

*  SATELLITES  TO  BE  SERVICED  IN  THE  ORBIT  DURING  THAT  MISSION.* 

*  8,  THE  MAXIMUM  8  OF  WAITING  ORBITS  TO  BE  MODELED.  THERE  IS  * 

*  ALSO  AN  OPTION  THAT  CALCULATES  PROPULSION  FUEL  MASS  USED.  * 


*  IF  IT  IS  DESIRED  TO  CALCULATE  PROPULSION  FUEL  MASS  USED. 

*  THE  USER  MUST  SPECIFY  THE  STRUCTURAL  MASS  OF  THE  OSV.  THE 

*  SPECIFIC  IMPULSE  OF  THE  OSV  FUEL.  AND  THE  AVERAGE  MASS 

*  OFFLOAD  AT  EACH  SATELLITE. 
*************************************************************** 

*  OUTPUT  OF  THIS  PROGRAM  IS  IN  THE  FORM  OF  TWO  CHARTS  WRITTEN* 

*  TO  TWO  FILES,  datal  AND  data2.  THE  CHART  IN  datal  GIVES  THE* 

*  FOLLOWING  INFORMATION  FOR  THE  8  OF  SATELLITES  SELECTED  TO  * 

*  BE  SERVICED  AND  THE  RANGE  OF  WAITING  ORBITS:  SEMI -MAJOR  * 

*  AXIS  OF  THE  WAITING  ORBIT.  PERIOD  OF  THE  WAITING  ORBIT.  * 

*  ALTITUDE  OF  WAITING  ORBIT  AT  PERIGEE  IN  KM  ABOVE  SURFACE  * 

*  OF  EARTH.  DELTA  V  OF  EACH  WAITING  ORBIT.  AND  DELTA  V  OF  * 

*  THE  RESUPPLY  ORBIT.  THE  CHART  IN  data2  SHOWS  THE  FOLLOWING  * 

*  INFORMATION  FOR  THE  8  OF  SATELLITES  SELECTED  TO  BE  SERVICED* 

*  AND  THE  RANGE  OF  WAITING  ORBITS:  TOTAL  DELTA  V  NEEDED  FOR  * 

*  THE  MISSION  SPECIFIED.  TIME  OF  FLIGHT  FOR  THAT  MISSION  AND  * 

*  PROPULSION  FUEL  MASS  BURNED  IN  ACCOMPLISHING  THE  MISSION  * 

*  (IF  FUEL  MASS  OPTION  IS  SELECTED).  * 

*************************************************************** 
*********  ************* 
*********  OUTPUT  WRITTEN  TO  'datal’  and  ’data2’  ************* 
*********  ************* 
*************************************************************** 
*********************  *************************** 
*********************  ASSUMPTIONS  *************************** 
*********************  *************************** 
*************************************************************** 

*  l)  SATELLITES  TO  BE  SERVICED  ARE  IN  CIRCULAR  ORBITS  AT  AN  * 
l  NCL  i  iv/ 

MODEL. 
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*************************************************************** 
C  *  2)  THE  NUMBER  OF  SATELLITES  IN  EACH  PLANAR  ORBIT  ARE  * 

C  *  EQUALLY  SPACED  ABOUT  THAT  ORBIT.  THIS  IS  A  CONSERVATIVE  * 
C  *  ESTIMATE  -  IF  SATELLITES  WERE  ACTUALLY  CLOSER.  DELTA  V 
C  *  REQUIRED  TO  SERVICE  EACH  SATELLITE  WOULD  BE  MUCH  LESS  THAN  * 
C  *  DELTA  V  NEEDED  IF  EQUALLY  SPACED.  CONSEQUENTLY  A  PLATFORM  * 
C  »  OF  SATELLITES  SHOULD  BE  TREATED  AS  A  SINGLE  SATELLITE  FOR  * 
C  *  ITS  DELTA  V  REQUIREMENTS.  * 

C  *************************************************************** 
C  *  3)  ALL  ORBIT  TRANSFERS  ARE  MODELED  AS  HOHMANN  TRANSFERS.  * 
C  *  WITH  TWO  IMPULSIVE  BURNS.  IF  A  PLANE  INCLINATION  CHANGE  * 
C  *  IS  INVOLVED.  THE  TRANSFER  IS  MODELED  AS  A  HOHMANN  TRANSFER  * 
C  *  WITH  PUNE  CHANGE  (PART  OF  CHANGE  AT  BURN  1  &  REST  AT  BURN  * 
C  *  2).  WAITING  ORBIT  ELLIPSES  (TRANSFER  ORBIT  OF  OSV  WAITING  * 
C  *  FOR  NEXT  SAT)  WILL  BE  ASSUMED  TO  BE  INSIDE  THE  CIRCULAR  * 
C  *  SERVICE  ORBIT.  TO  ALLOW  A  MINIMUM  WAITING  TIME  TO  BE  USED.  * 
C  *  HOWEVER.  THE  PERIGEE  OF  THE  WAITING  ORBIT  MUST  BE  GREATER  * 
C  *  THAN  THE  RADIUS  OF  THE  EARTH  +  IOO  NM  TO  NEGLECT  MAJOR  * 
C  *  ATMOSPHERIC  DRAG  EFFECTS  (AND  TO  AVOID  RUNNING  INTO  THE  * 
C  *  EARTH).  R( PERIGEE. WO)  >  REARTH  -*•  100  NM  * 

C  *  * 

C  *  THE  CHART  IN  data I  FILE  SHOWS  THE  ALTITUDE  ABOVE  THE  EARTHS* 
C  -  SURFACE  AT  PERIGEE.  SO  THAT  YOU  CAN  SELECT  THE  NUMBER  OF  * 
C  *  WAITING  ORBITS  NECESSARY  TO  KEEP  PERIGEE  ABOVE  100  NM.  * 
C  *************************************************************** 
C  *  4)  SERVICING  IS  ASSUMED  TO  BE  REGUURLY  SCHEDULED  PRE-  * 
C  *  VENTATIVE  MAINTENANCE.  CONSEQUENTLY  EACH  SATELLITE  WILL  * 
C  *  BE  SERVICED  IN  CONSECUTIVE  ORDER  WITH  AN  AVERAGE  AMOUNT  * 
C  *  OF  MASS  OFFLOADED  AT  EACH  SATELLITE.  * 

C  *************************************************************** 
C  *************************************************************** 

C  *  DELTA  V  REQUIREMENTS  ARE  DETERMINED  BY  THE  ORBITAL  * 

C  *  MECHANICS  INVOLVED.  ONCE  DELTA  V  IS  OBTAINED.  FUEL  MASS  * 
C  *  REQUIREMENTS  MAY  BE  OBTAINED  USING  THE  ROCKET  EQUATION:  * 

C  *  DELTA  V  =  lsp*G*ln(Mass  initial/Mass  final)  * 

C  where  Mass  initial  =  Mass  structure  *  Mass  fuel  +  Mass  payload 
C  *  * 

C  *  and  Mass  final  =  Mass  structure  Mass  payload  * 

C  *  * 

C  *  DELTA  V  TOTAL/MJSSION  =  (DELTA  V  required  from  rendezvous  * 

C  *  with  supplies  to  service  orbit  and  * 

C  *  return)  +  (DELTA  V  required  from  * 

C  *  service  orbit  to  waiting  orbit  and  * 

C  *  return )*( number  of  satellites  ser-  * 

C  *  viced  -  1 )  » 

C  *  * 

C  *  * 

C  *  ONE  MISSION  IS  DEFINED  AS  TIME  FROM  BEING  SUPPLIED  WITH  * 
C  *  MASS  TO  RETURN  FOR  RESUPPLY  AFTER  SERVICING  SATELLITES.  * 
C  *  * 

C  *  PARKING  ORBIT  IS  DEFINED  AS  THE  ORBIT  WHERE  RENDEZVOUS  * 
C  *  WITH  RESUPPLIES  IS  ACCOMPLISHED.  * 

C  * 
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c  * 

WAITING  ORBIT  IS 

DEFINED  AS  THE  ORBIT  USED  BY 

THE 

osv 

* 

c  * 

WHEN  WAITING  FOR 

THE  NEXT  SATELLITE  TO  "CATCH 

UP" 

TO 

■* 

c  * 

BE  SERVICED. 

* 

C  **************************************************************** 
c  *************************************************************** 

c  *********#*)HHH()Ht#**lt)HHH(##)(*##iH(*iHHHH<1HHHHt*#*******i(****lH(*#* 

C  *  THE  FOLLOWING  RULES  ARE  USED  FOR  DOING  INCLINATION  CHANGES:  * 
C  *  (HOHMANN  TRANSFER  WITH  PLANE  CHANGE  WITH  2  IMPULSIVE  BURNS)  * 


C  *  * 

C  *  1)  IF  THE  INCLINATION  CHANGE  IS  LESS  THAN  OR  EQUAL  TO  ONE  * 

C  *  DEGREE.  THE  ENTIRE  CHANGE  IS  ACCOMPLISHED  DURING  THE  * 

C  *  SECOND  BURN.  * 

C  *  * 

C  *  2)  IF  THE  INCLINATION  CHANGE  IS  LESS  THAN  OR  EQUAL  TO  FIVE  * 

C  *  DEGREES.  ONE  DEGREE  OF  PLANE  CHANGE  IS  DONE  DURING  THE  * 

C  *  FIRST  BURN,  AND  THE  REMAINDER  DURING  THE  SECOND  BURN.  * 

C  *  * 

C  *  3)  FOR  PLANE  CHANGES  GREATER  THAN  FIVE  DEGREES.  THREE  * 

C  *  DEGREES  OF  PLANE  CHANGE  IS  DONE  DURING  THE  FIRST  BURN.  * 

C  *  AND  THE  REMAINDER  DURING  THE  SECOND  BURN.  * 

C  *  * 

C  *  THESE  RULES  WERE  ARRIVED  AT  AFTER  CAREFUL  COMPUTER  * 

C  *  MODELING  AND  ANALYSIS  AND  WERE  FOUND  TO  VARY  FROM  THE 

C  *  OPTIMAL  BY  LESS  THAN  1%.  * 


C  *************************************************************** 
C  *************************************************************** 
C  *************************************************************** 
double  precision  two. a. vso.vwoa.de I vwo . vtra , vpo , vtrp , 

-t-tdelv  .mi  , fm.mf  .delvp.delva.delvrs 
real  mu.rearth.altsat.rsat.pi .alt.raddif , 
i-altpo.  i . n  .a  1  ti t . wmax . smax , rad , radl  ,rad3, 

+osvms. isp.mpld.gc.sat. inc, il . i2 
integer  ans 

open  (unit=3.f i le= 'datal ’ ) 
open  (unit=4.file=’data2’ ) 
mu = 3 . 9860 1 2E^05 

c  *********  Mil  is  GRAVITATIONAL  PARAMETER  OF  EARTH  IN  KM  ******** 
rearth=6378. 145 

c  «****#«*  REARTH  IS  RADIUS  OF  EARTH  IM  KM  ********************* 
write  (*.*)  ’what  is  altitude  in  km  of  satellite  class'. 

’  you  want  to  transfer  to?  ’ 
read  (*,*)  altsat 

write  (*,*)  ’what  is  the  inclination  in  degrees  of  the', 
-•-’satellite  class  you  want  to  transfer  to?’ 
read  (*,*)  i2 

c  *********  ALTSAT  IS  ALTITUDE  OF  SATELLITE  CLASS  IN  KM  ********* 
write  ( 3,*j  'PROGRAM  NAME:  vorbcs.f’ 
write  (3,*)  ’EXECUTABLE  FILE:  vorbcs’ 
write  (3.*)  ' ’ 

write  (3.*) ’THIS  IS  FILE  datal,  ADDITIONAL  INFORMATION  MAY', 
+’BE  FOUND  IN  FILE  data2.’ 

WI1  tc  (3  ,  S  ) 

write  (4.*)  ’PROGRAM  NAME:  vorbcs.f’ 
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write  (.4.*)  EXECUTABLE  FILE.  vorbcs’ 
write  (4.*;  ' ’ 

write  !4  ,*)  'THIS  IS  FILE  data2.  ADD  I T I ONAL  INFORMATION  ’  . 

-  ’MAY  BE  FOUND  IN  FILE  data!.’ 
write  (4 .*)  ’ ' 

write  (3.*)  ’DOUBLE  PRECISION  MODEL’ 
write  (4.*)  DOUBLE  PRECISION  MODEL’ 

write  !*.*)  ’ALTITUDE  OF  SATELLITE  CLASS  IS  ’ .altsat.  KM’ . 

-  ’  AND  INCLINATION  OF  CLASS  IS’ . i2. ’DEGREES’ 

write  (3.*)  ’ALTITUDE  OF  SATELLITE  CLASS  IS  ’ .altsat. ’KM’ . 
-•-’AND  INCLINATION  OF  CLASS  IS’  . i2.  ’ DEGREES’ 
write  (4.*)  ’ALTITUDE  OF  SATELLITE  CLASS  IS  ’ .altsat. ’KM’ . 
-AND  INCLINATION  OF  CLASS  IS’ .i2. ’DEGREES’ 
rsat=a l tsat-rearth 
write!*,*)  ’ ’ 

write  !*.*)’what  is  altitude  in  km  of  parking  orbit  for’, 
-’rendezvous  with  resupplies?’. 

-  ’  ( 1  OOnm=  185.2.1 50nm= 277.8. 200nm=  370.4)’ 

c*altpo  is  altitude  of  parking  orbit  for  rz  with  resupplies  in  km 

read  ( * .  * )  a  1  tpo 

write!* .*) ’what  is  the  inclination  in  degrees  of  the’, 
-•-’parking  orbit?’ 
read!*.*)  il 

write  !*.*)  ’ALTITUDE  OF  PARKING  ORBIT  IS’.altpo.'  KM’, 
-’AND  INCLINATION  IS’  .il .  DEGREES’ 

write  !3.*)  ALTITUDE  OF  PARKING  ORBIT  IS’.aitpo.’  KM’. 
-AND  INCLINATION  IS ’. i 1 , ’DEGREES’ 

write  (4.*)  ALTITUDE  OF  PARKING  ORBIT  IS’.altpo.’  KM’. 
-AND  INCLINATION  IS’ . il .  DEGREES’ 
a  1 t=a 1 tpo-rearth 
inc=abs! 12— i 1 ) 
pi =3. 1415927 
write  !*.*)’’ 

write  !*,*)  ’enter  the  8  of  satellites  in  the  orbit’. 

-’to  which  you  are  transferring’ 
read  ! *  .  * )  smax 

write  !*.*)  '8  OF  SATELLITES  IN  ORBIT  IS’. smax 
write!*,*) ’ ’ 

write  (*.*)  ’enter  max  8  of  orbits  you  want  modeled.’ 

-’for  waiting  orbit’ 
write  !*.*)’’ 

write  !*.*)' NOTE :  8  OF  WAITING  ORBITS  IS  THE  NUMBER'. 

-OF  REVOLUTIONS  OSV  MAKES  IN  TRAVELING  FROM  SATELLITE  TO’. 
-’SATELLITE.  FOR  MORE  WAITING  ORBITS.  LESS  DELTA  V  AND’. 
-’PROPULSION  FUEL  IS  USED.  BUT  MORE  TIME  IS  REQUIRED.’ 
read  !*,*)  wmax 

write  !*.*)  ’MAX  8  OF  ORBITS  MODELED  IN  WAITING  ORBIT’. 
-’WILL  BE ’ . wmax 
write!*,*)  ’ ’ 

write  !*.*)’how  many  satellites  do  you  want  to  service’ . 
-on  one  mission?’ 
read!*,*)  sat 
write  !*.*)  ’ ’ 
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write!* .*) 'THIS  PROGRAM  WILL  CALCULATE  TOTAL  DELTA  V.'. 
-’TIME  OF  FLIGHT.  AND  (IF  DESIRED)  PROPULSION  FUEL  MASS’. 
-’NECESSARY  FOR  A  MISSION  SERVICING’ .sat.  SATELLITES  OUT’. 
-’OF  A  TOTAL  CONSTELLATION  OF ’. smax .’ SATELL I TES ’ 
write!* ,*)  ’  ’ 
write  (3.*)  ’  ’ 

write (3.*) 'THIS  PROGRAM  WILL  CALCULATE  TOTAL  DELTA  V.’. 
-’TIME  OF  FLIGHT.  AND  (IF  DESIRED)  PROPULSION  FUEL  MASS’. 
-’NECESSARY  FOR  A  MISSION  SERVICING’ .sat. ’SATELLITES  OUT’. 
-OF  A  TOTAL  CONSTELLATION  OF’ .smax.  SATELLITES’ 
write(3.*)  ’’ 
write  (4 , *)  ’ ’ 

write (4.*)  THIS  PROGRAM  WILL  CALCULATE  TOTAL  DELTA  V.’. 
-’TIME  OF  FLIGHT.  AND  (IF  DESIRED)  PROPULSION  FUEL  MASS’. 
-’NECESSARY  FOR  A  MISSION  SERVICING’ .sat. 'SATELLITES  OUT’. 
-’OF  A  TOTAL  CONSTELLATION  OF ’.smax. ’SATELLITES’ 
write(4 . *)  ’ ’ 

write!*.*) ’do  you  want  to  calculate  osv  propulsion  fuel’, 
-’burned?  answer  with  number:  yes=I  or  no=0’ 
read  (*.*)  ans 
if  (ans.eq.O)  go  to  400 
write(*.*)  ’ ’ 

write  (*.*)  ’What  is  the  empty  structural  mass  of  your’. 
-OSV  in  kg?’ 
read  ( * . * )  osvms 

write  (*.*)  'OSV  STRUCTURAL  MASS  IS’ .osvms.  KG’ 
write  (3.*)  ’OSV  STRUCTURAL  MASS  IS’  ,->svms .  KG’ 
write  (4.*)  OSV  STRUCTURAL  MASS  IS ’ .osvms . ’ KG ’ 
write(*.») ’ ’ 

write  (*.*)  ’What  is  the  specific  impulse  of  the  OSV’, 
-’fuel  in  sec?’ 
read  ( * . * )  i sp 

write  (*,*)  ’OSV  FUEL  Isp  =’,isp,’SEC’ 
write  (3.*)  ’OSV  FUEL  Isp  =',isp.'SEC’ 
write  (4.*)  ’OSV  FUEL  Isp  =’.isp,’SEC’ 
write(*.*)  ’ * 

write  (*.*)  ’What  is  the  average  mass  offload  at  a’, 
-’satellite  in  kg?’ 
read  ( * , * )  mp 1 d 

write  (*.*)  ’AVE  PAYLOAD  OFFLOAD  AT  EACH  SATELLITE  =’. 

-mp 1 d . ’ KG ’ 

write  (3.*)  ’AVE  PAYLOAD  OFFLOAD  AT  EACH  SATELLITE  =’. 

-mpld.  KG’ 

write  (4.*)  ’AVE  PAYLOAD  OFFLOAD  AT  EACH  SATELLITE  =’. 
-mpld. ’KG’ 

write(*.*) ’’PLEASE  STANDBY  -  I’M  BUSY  COUNTING! !" 

C  ****  gc  is  acceleration  of  gravity  *********** 
gc=9.8G7E-03 

write(3.*) ’FUEL  MASS  IS  TOTAL  OSV  FUEL  USED  FOR  PROPULSION ’ 
-IN  KG’ 

write(4.») ’FUEL  MASS  IS  TOTAL  OSV  FUEL  USED  FOR  PRUPULSION’ 
-  IN  KG 

400  write(3.*)  PERIOD  OF  THE  WAITING  ORBIT  IS  IN  HOURS’ 
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'CROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  bureau  of  standards  1%3  a 


writeU.*)  PERIOD  OF  THE  WAITING  ORBIT  IS  IN  HOURS’ 
write  (3.*) ’ ’ 

write(3.*) ’THE  a  OF  WAITING  ORBITS  INDICATES  THE  NUMBER  OF’ 
-REVOLUTIONS  REQUIRED  TO  TRANSIT  FROM  ONE  SATELLITE  TO  THE’ 
- ’ NEXT  NOTE  THAT  THE  MORE  WAITING  ORBITS  USED  BETWEEN’. 
-’SATELLITES.  THE  LESS  DELTA  V  AND  PROPULSION  FUEL  NEEDED.’ 
write(3.*) ’ ’ 
write  (4,**)  ’  ’ 

write(4.*) ’THE  a  OF  WAITING  ORBITS  INDICATES  THE  NUMBER  OF’ 
-'REVOLUTIONS  REQUIRED  TO  TRANSIT  FROM  ONE  SATELLITE  TO  THE’ 
-’NEXT  NOTE  THAT  THE  MORE  WAITING  ORBITS  USED  BETWEEN’. 
-’SATELLITES,  THE  LESS  DELTA  V  AND  PROPULSION  FUEL  NEEDED.’ 
write(4.*) ’ ’ 
write!*.*) ' ’ 

write!*.*)  ’’PLEASE  STANDBY  -  I’M  BUSY  COUNTING!!” 
write  (3.*)  ’DELTA  V  IS  IN  KM/HR’ 
write  (3.*)  ’ ’ 

write  (4.*)  DELTA  V  IS  IN  KM/HR’ 
write  (4.*)  ’ ’ 

wr  i  te  (  3 .  * )  ’  *hhwhhmhhhmwhhhhhhhhmhhhmhhhhhmhhmhhhhhhhhhhh*#  ’ 
write  ( 4 , * )  ’  **********************************************  ’ 
write  (3,*) ’ ’ 
write  (4,*) ’ ’ 
write  (3.10) 

10  format! ’a  OF  SATS’,4X.’»  OF  SATS’.4X.’«  OF  WAITING  ORBITS’, 
+6X,  ’PERIOD. WO’  ,10X.  'A, WO’  ,10X.  *W0  ALT  ABV\8X, ’DELTA  V.WO’. 
-5X. ’DELTA  V. RESUP’) 

write  (3.11) 

11  format! lx, 'SERVICED' .4X, ’ IN  PLANE’ ,33X. 'HRS’ . 14x. ’ KM’ . 13X . 
-’EARTH.  KM’ ,10X,’ KM/HR’ , 11X,  KM/HR’ ) 

write  (3.*)  ’ ’ 

IF  (ans.eq.O)  go  to  20 
write  (4.25) 

25  format  (’«  OF  SATS’.4x,’#  OF  SATS’,4X,’a  OF  WAITING  ORBITS’ 
-4x,  TOTAL  DELTA  V’,4X.’TIME  OF  FLIGHT’ ,4X. ’FUEL  MASS’) 

write  (4,12) 

12  format( lx.  SERVICED’ .4X. 'IN  PLANE' .31X, ’KM/HR’ . 13X. ’HRS’ , 
-12x,  KG’) 

write  (3,*)  ’ ’ 
go  to  30 
20  write  (4.26) 

26  format  (’«  OF  SATS’,4x.’«  OF  SATS'.4X.’a  OF  WAITING  ORBITS’ 
+4x. 'TOTAL  DELTA  V’.4X,’TIME  OF  FLIGHT’) 

write  (4.13) 

13  formatdx. ’SERVICED’ .4X, ’IN  PLANE’ ,31X. ’KM/HR’ .13X. ’HRS’ ) 
write  (3.*)  ’ ’ 

30  do  100  nssmax,smax 

C  *********  N  IS  NUMBER  OF  SATELLITES  PER  ORBIT  ************ 
do  200  i = l , wmax 

C  *********  I  is  NUMBER  OF  ORBITS  IN  WAITING  ORBITS  ************* 
two=( (n*i-l )/ (n*i ) )*(rsat*»(3. /2. ) )*2.*pi/sqrt(mu) 

C  *********  TWO  IS  PERIOD  OF  WAITING  ORBIT  ********************** 
ttro=2.*pi*sqrt((((alt+rsat)/2. )**3. )/mu) 


c  *********  TTRO  IS  PERIOD  OE  TRANSFER  ORBIT  (PO  TO  SO)  ********* 
tof =ttro-i-two*( sat-1 . ) 

C  *********  TOF  IS  TIME  OF  FLIGHT  FOR  ONE  COMPLETE  MISSION  ****** 
a=(two*sqrt(mu)/(2.*pi))**(2./3. ) 

C  *********  a  IS  SEMI -MAJOR  AXIS  OF  WAITING  ORBIT  *************** 
altit=2.*a-rsat-rearth 

C  *********  ALT1T  IS  ALTITUDE  ABOVE  EARTH  AT  PERIGEE  OF  WO  ******* 
vso= sqrt ( mu / r sat ) 

C  *********  VSO  IS  SPEED  IN  SERVICE  ORBIT  ************************ 
vwoa=sqrt(2.*(rou/rsat-mu/(2.*a) )) 

C  *********  VWOA  IS  SPEED  IN  WAITING  ORBIT  AT  APOGEE  ************ 
de 1 vwo=2 . *dabs ( vso- vwoa ) 

C  **  DELVWO  IS  DELTA  V  FOR  TRAVELING  BETWEEN  SATS  (WAITING  ORBIT) 
vtra=sqrt(2.*abs(mu/rsat-mu/ (alt+rsat) ) ) 

C  ****  VTRA  IS  SPEED  IN  TRANSFER  ORBIT  AT  APOGEE  FROM  SO  TO  PO  ** 
vpo=sqrt(mu/ait) 

C  ********  VPO  IS  VELOCITY  OF  CIRCULAR  PARKING  ORBIT  ********* 
vtrp=sqrt(2.*abs(mu/alt-mu/(alt+rsat) ) ) 

C  *  VTRP  IS  SPEED  OF  OF  TRANSFER  ORBIT  AT  PERIGEE  FROM  PO  TO  SO  * 

C 

C  *************************************************************** 

C  *  THE  FOLLOWING  RULES  ARE  USED  FOR  DOING  INCLINATION  CHANGES:  * 


C  *  * 
C  *  1)  IF  THE  INCLINATION  CHANGE  PS  LESS  THAN  OR  EQUAL  TO  ONE  * 
C  *  DEGREE.  THE  ENTIRE  CHANGE  IS  ACCOMPLISHED  DURING  THE  * 
C  *  SECOND  BURN.  * 
C  *  * 
C  *  2)  IF  THE  INCLINATION  CHANGE  IS  LESS  THAN  OR  EQUAL  TO  FIVE  * 
C  *  DEGREES.  ONE  DEGREE  OF  PLANE  CHANGE  IS  DONE  DURING  THE  * 
C  *  FIRST  BURN.  AND  THE  REMAINDER  DURING  THE  SECOND  BURN. 

C  *  * 
C  »  3)  FOR  PLANE  CHANGES  GREATER  THAN  FIVE  DEGREES.  THREE  * 
C  *  DEGREES  OF  PLANE  CHANGE  IS  DONE  DURING  THE  FIRST  BURN.  * 
C  *  AND  THE  REMAINDER  DURING  THE  SECOND  BURN.  * 
C  *  * 
C  *  THESE  RULES  WERE  ARRIVED  AT  AFTER  CAREFUL  COMPUTER  * 
C  *  MODELING  AND  ANALYSIS  AND  WERE  FOUND  TO  VARY  FROM  THE  * 
C  *  OPTIMAL  BY  LESS  THAN  1%.  * 
C  *************************************************************** 

c  *  * 

C  *  RAD  IS  THE  PLANE  INCLINATION  DIFFERENCE  IN  RADIANS  * 


rad=inc*pi/ 180. 

C  *  RADI  IS  ONE  DEGREE  IN  RADIANS  **** 
radl=pi/180. 

C  *  RAD3  IS  THREE  DEGREES  IN  RADIANS  **** 
rad3=pi/60. 

IF  (inc.le.1.0)  THEN 

de I vp=dsqr t ( ( vpo**2 .  0 ) + ( v tr p**2 . 0 ) -2 . 0* vpo*vt rp ) 
delva=dsqrt((vso**2.0)+(vtra**2.0)-2.0*vso*vtra*cos(rad) ) 
ELSE  IF  (inc.Ie.5.0)  THEN 

delvp=dsqrt( (vpo**2.0)+(vtrp**2.0)-2.0*vpo*vtrp*cos(radl ) ) 
raddif=rad-radl 

delva=dsqrt( (vso**2.0)+(vtra**2.0)-2.0*vso*vtra*cos(raddif ) ) 
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ELSE 

de i vp=dsqrt ( ( vpo**2 . 0 )+ ( vtrp**2 . 0 )-2 . U*vpo*vtrp*cos ( rad3 ) ) 
raddif=rad-rad3 

de  1  va=dsqrt  ( ( vso**2 .  0 ) + (  vtra**2 . 0 )  -2 .  O*vso*vtra*cos  (  radd  i  f  > ) 
END  IE 

delvrs=2.*(delvp-i-deiva) 

C  ******  DELVRS  IS  DELTA  V  FOR  RESUPPLY  ******** 
tdelv=delvwo*(sat~l )+delvrs 
C  **  TDELV  IS  TOTAL  DELTA  V  ** 
if  (ans.eq.O)  go  to  70 

C  **  mi  IS  MASS  INITIAL,  mf  IS  MASS  FINAL,  fm  IS  FUEL  MASS  ****** 
mi =osvms*dexp (delvrs/(2.*i sp*gc ) ) 
fm=mi-osvms 
mf=mi+mpld 

IF  (sat.eq.l.)  go  to  300 
do  300  j=l .sat- I . 
mi=mf*dexp(del vwo/ ( isp*gc) ) 
fm=fm-^mi-mf 
mf-mi^mpld 
300  continue 

mi=raf*dexp(delvrs/ (2.*isp*gc) > 
fm=fm+rai-mf 

C  ******  CHANGE  TIME  IN  ANSWERS  FROM  SEC  TO  HOURS  *********** 

70  two=two/3600. 

de 1 vwo=de 1 vwo*3600 . 
del vrs=del vrs*3600 . 
tde 1 v=tde l v*3600 . 
tof=tof /3600. 

write  (3.50)sat,n. i .two.a.altit.delvwo.delvrs 
50  format  (F5.0.8x.F5.0,l5x,F5.0.7x.5Dl7.9) 

IF  (ans.eq.O)  go  to  60 

write  (4.75)  sat.n. i .tdelv.tof ,fm 

75  format  (F5.0.8x .F5.0, 15x .F5.0.8x .3D17 .9) 
go  to  200 

60  write  (4,76)  sat.n. i .tdelv.tof 

76  format  (F5.0.8x.F5.0. l5x.F5.0.8x.2D17.9) 

200  continue 

100  continue 

write  (*.*)  "INFORMATION  IS  IN  FILE  data!’  and  data2” 
write  (3.*)  ’ ' 
write  (4.*)  ’ ’ 

write  (3.*) ’*********************************************’ 
write  (4 . *) ’ ********************************************* ’ 
end 
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PROGRAM  NAME:  vorbcs.f 
EXECUTABLE  FILE:  vorbcs 

THIS  IS  FILE  datal .  ADDITIONAL  INFORMATION  MAY  BE  FOUND  IN  FILE  data2. 

DOUBLE  PRECISION  MODEL 

ALTITUDE  OF  SATELLITE  CLASS  IS  800.000  KM  AND  INCLINATION  OF  CLASS  IS  0. 
DEGREES 

ALTITUDE  OF  PARKING  ORBIT  IS  185.200  KM  AND  INCLINATION  IS  0.  DEGREES 

THIS  PROGRAM  WILL  CALCULATE  TOTAL  DELTA  V, 

TIME  OF  FLIGHT.  AND  (IF  DESIRED)  PROPULSION  FUEL  MASS 
NECESSARY  FOR  A  MISSION  SERVICING  5.00000  SATELLITES  OUT 
Or  A  TOTAL  CONSTELLATION  OF  144.000  SATELLITES 

0S\  STRUCTURAL  MASS  IS  2000.00  KG 
OSV  FUEL  Isp  =  450.000  SEC 

AVE  PAYLOAD  OFFLOAD  AT  EACH  SATELLITE  =  250.000  KG 

FUEL  MASS  IS  TOTAL  OSV  FUEL  USED  FOR  PROPULSION  IN  KG 
PERIOD  OF  THE  WAITING  ORBIT  IS  IN  HOURS 

THE  ?  OF  WAITING  ORBITS  INDICATES  THE  NUMBER  OF 
RESOLUTIONS  REQUIRED  TO  TRANSIT  FROM  ONE  SATELLITE  TO  THE 
NEXT  NOTE  THAT  THE  MORE  WAITING  ORBITS  USED  BETWEEN 
SATELLITES.  THE  LESS  DELTA  V  AND  PROPULSION  FUEL  NEEDED. 

DELTA  V  IS  IN  KM/HR 


s  OF  SATS 
SERVICED 


8  OF  SATS 
IN  PLANE 


a  OK  WAITING  ORBITS 


PERIOD. WO 
HRS 


A.  WO 
KM 


144. 


1-  0. 166955200d+01  0.7l4187395d->04 


*  *  *  **  *»»**«#*««***#*****»»*******#*.# 


V 
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LO  ALT  ABV 

earth,  km 

0 .733457B86d+U3 


UO  ALT  ABV  DELTA  V.WO  DELTA  V.KESUP 

EARTH,  KM  KM/HR  KM/ HR 

733457885d-^03  0. 125065556d+03  0.245542030d+04 


PROGRAM  NAME:  vorbcs.f 
EXECUTABLE  FILE:  vorbcs 

THIS  IS  FILE  data2.  ADDITIONAL  INFORMATION  MAY  BE  FOUND  IN  FILE  datal . 

DOUBLE  PRECISION  MODEL 

ALTITUDE  OF  SATELLITE  CLASS  IS  800.000  KM  AND  INCLINATION  OF  CLASS  IS  0. 
DEGREES 

ALTITUDE  OF  PARKING  ORBIT  IS  185.200  KM  AND  INCLINATION  IS  0.  DEGREES 

THIS  PROGRAM  WILL  CALCULATE  TOTAL  DELTA  V, 

TIME  OF  FLIGHT,  AND  (IF  DESIRED)  PROPULSION  FUEL  MASS 
NECESSARY  FOR  A  MISSION  SERVICING  5.00000  SATELLITES  OUT 
OF  A  TOTAL  CONSTELLATION  OF  144.000  SATELLITES 

OSV  STRUCTURAL  MASS  IS  2000.00  KG 
OSV  FUEL  Isp  =  450.000  SEC 

AYE  PAYLOAD  OFFLOAD  AT  EACH  SATELLITE  =  250.000  KG 

FUEL  MASS  IS  TOTAL  OSV  FUEL  USED  FOR  PROPULSION  IN  KG 
PERIOD  OF  THE  WAITING  ORBIT  IS  IN  HOURS 

THE  a  OF  WAITING  ORBITS  INDICATES  THE  NUMBER  OF 
REVOLUTIONS  REQUIRED  TO  TRANSIT  FROM  ONE  SATELLITE  TO  THE 
NEXT  NOTE  THAT  THE  MORE  WAITING  ORBITS  USED  BETWEEN 
SATELLITES.  THE  LESS  DELTA  V  AND  PROPULSION  FUEL  NEEDED. 

DELTA  V  IS  IN  KM/ HR 


5  OF  SATS  a  OF  SATS  s  OF  WAITING  ORBITS  TOTAL  DELTA  V  TIME  OF  FLIGHT  FUEL  MASS 
SERVICED  IN  PLANE  KM.  HR  HRS  KG 

144-  1-  0 . 2955b8253d+04  0 .825260353d+01  0.5808808b7d^0 


*»»*mU«*»K)H(*)(»***»****«*|Ht»»|f»»*»*»»»*l(***» 


U- 


H-10 


Bi bLiogcaeby 


"Aerospace  Firms  Offer  OTV  Designs."  Space  World.  A-2-218: 
24-30-1-  (February  1982). 

Aldridge.  E.  C.  Jr.  Memorandum  for  the  Vice  Chief  of  Staff. 
USAF.  Department  of  the  Air  Force.  Washington  DC.  13 
September  1984. 

Baker.  Jerome  M.  Orbit  Transfer  and  Rendezvous  Mannevers 
Between  Inclined  Circular  Orbits.  Contract  No.  AF 
04(695)-469.  Aerospace  Corporation.  El  Segundo  CA ,  June 
1965  (SSD-TR-65-7b) . 

Bate.  Roger  R.  et  a  1  .  Fundamentals  of  Astrodvnamics .  New 
York:  Dover  Publications.  Inc..  1971. 

Bolen.  Major  C.  T.  et  at.  Lunar  Supply  Shuttle  Preliminary 
Des i gn .  Graduate  Astronautics  Class  Project  in  MC  b.7b. 
School  of  Engineering.  Air  Force  Institute  of  Technology 
(AU),  Wr i ght-Patter son  AFB  OH.  May  1980 

Briggs.  Hugh  C.  and  Aaron  R.  DeWispelare.  I nteract  i  ve 
Analysis  and  Interpretive  Structural  Modelling  in  a 
Program  Planning  and  System  Design.  Unpublished 
Technical  Report  No.  AU-AF 1 T-EN-TR-X .  School  of 
Engineering.  Air  Force  Institute  of  Technology,  Wrignt 
Patterson  AFB  OH.  1985a. 

Carlton.  Robert.  Senior  Staff  Engineer.  Aerospace 

Corporation.  Personal  interview.  NASA  Johnson  Space 
Center.  Houston  TX .  23  August  1985. 

Chankong.  Vira  and  Yacov  Y.  Haimes.  Mu  1 1  iob.iect  i  ve  Decision 
Making:  Theory  and  Methodology.  New  York:  Elsevier 
Science  Publishing  Company.  Inc..  1983 

Chestnut,  Harold.  Systems  Engineering  Tocls.  New  York- 
Wiley  and  Sons,  Inc..  1965. 

Clark.  Dennis  R.  and  Aaron  R.  DeWispelare.  "Application  of 
Multiple  Objective  Optimization  Theory  to  a  Preliminary 
Missile  Design  Study.”  Design  Studies.  6-2:  (April 
1985). 

- .  User’s  Manual  for  Solving  Nonlinear  Multiple 

Objective  Problems  Using  Optimization  Routine  SUMT 
(Program  Proces).  Wright  Patterson  AFB  OH:  Air  Force 
Institute  of  Technology.  1983. 


BIB-1 


Covault.  Craig.  ’’Aeronautics  Policy  Stresses  Transport. 
Transatmospher i c  Work."  Aviation  Week  and  Space 
Technology .  122  ■  lb-20---  (April  1985). 

Crabtree.  Colonel  W.H..  Deputy  Commander  for  Space  Systems. 
USAF  Space  Division.  Personal  interview.  Air  Force 
Institute  of  Technology,  WPAFB  OH.  20  August  1985. 

Crawford.  Gordon  and  Cindy  Williams.  The  Ano lysis  of 
Subjective  Judgement  Matrices.  Project  AIR  FORCE 
Report,  Rand  Corporation.  Santa  Monica  CA .  May  1985 
(R-2572-1-AF) . 

DeNeuf v i 1 1 e  .  Richard  and  Joseph  H.  Stafford.  Systems 

Analysis  for  Engineers  and  Managers.  New  York:  McGraw- 
Hill  Book  Company.  Inc..  1971. 

Department  of  the  Air  Force.  Space  Handbook.  AU~18. 

Alabama:  Maxwell  AFB .  August  1985. 

Headquarters  Space  Division  (AFSC).  AFSC  Field 
Activity  Management  Policy.  Space  Servicing.  SD 
Regulation  540-8.  Los  Angeles:  HQ  SD.  25  October  1984a. 

USAF  Shuttle  Assessment  Report.  Report  SD-CF-0192. 
Los  Angeles:  USAF  Space  Division.  October  1984b. 

- •  United  States  Air  Eorce  Spacecraft  Maintenance  Policy 

Rev iew .  Volume  I.  Study  Report.  HQ  USAF.  Washington  DC, 
June  1984c. 

- .  United  States  Air  Force  Spacecraft  Maintenance  Policy 

Review .  Volume  II.  Supporting  Documentation.  HQ  USAF, 
Washington  DC,  June  I984d. 

- .  Space  Planners  Guide.  Washington:  Goverment  Printing 

Office.  1965. 

DeWispelare.  Aaron  A.  Vector  Optimization  Techniques. 

Technical  Report  No.  AU -AF 1 T-EN-TR-83-5  .  School  of  . 
Engineering.  Air  Force  Institute  of  Technology.  Wright 
Paterson  AFB  OH.  1984. 

- .  The  Extended  Multi-Attribute  Decision  Analysis  Mode  1 

( EMADAM ) .  Technical  Report  No.  AU-AF 1 T-EN-TR-83-3  . 

School  of  Engineering.  Air  Force  Institute  of 
Technology.  Wright  Paterson  AFB  OH.  1983. 

DeWispelare.  Aaron  R.  and  Wayne  A.  Stimpson.  User's  Manual 
for  program  MADAM  (The  Multiple  Attribute  Decision 
Ana  lysis  Mode  1 ) .  School  of  Engineering.  Air  Force 
Institute  of  Technology.  Wr ight-Patterson  AFB,  OH.  1983. 


BI  P-2 


Dooling.  Dave.  "Waiting  for  the  Repairman,-’  Space  World,  S__ 
2~218  ■  4-35-r  (February  1982). 

Eklund.  Philip.  "Space-Borne  Projectile."  L  5  News  ,  9-5  ■  iT 
(May  1984). 

Feldman.  Major  James  K.  and  Major  William  F.  Rowell.  Lecture 
materials  distributed  in  SM  b . 33 ,  Systems  Effectiveness 
Trade-off  Studies.  School  of  Engineering,  Air  Force 
Institute  of  Technology.  Wr ight-Patterson  AFB.  OH. 
January  1985. 

Fiaco.  Anthony  V.  and  Garth  Me  Cormick.  toon- L i nea r 

Programming  Sequential  Unconstrained  Minimization 
Techn i ques .  Mew  York:  John  Wiley  and  Sons.  1968. 

Fong.  Franklin  K.  et  a  1 .  Unmanned  Spacecraft  Cost  Model. 
Technical  Report  Mo.  SD-TR-81-45.  Directorate  of  Cost 
Analysis.  Space  Division.  Los  Angeles.  CA .  1981. 

Franklin.  J.N.  Matrix  Theory.  Englewood  Cliffs,  MJ : 
Prentice-Hall.  1968. 

Frisbee.  Robert  H.  ’’Propulsion  Systems-1,”  Space  Education. 

1  394-400  (9  May  1985).  - 

General  Dynamics /Conva i r  Division,  Presentation  Materials. 
Economic  Analysis  of  Satellite  Servicing.  Committee  on 
Spacecraft  Maintenance  and  Repair.  NASA  Headquarters. 
Washington  D.C..  June  1983. 

Goicoechea.  Ambrose,  et  a  1  .  Mu  1 1  iob.iect  i  ve  Decision  Analysis 
with  Engineering,  and  Business  Applications.  Mew  York: 
John  Wiley  and  Sons.  Inc..  1982. 

Green.  Col  Gaylord.  Director  of  Plans.  Personal  interview. 

HO  SD.  Los  Angeles  AFS  CA.  2  July  1985. 

Guy.  Walter.  "Partially  Closed  System  Will  Save  Big  Money.” 
Astronautics  &  Aeronautics.  March  1983. 

Haimes.  Yacov  Y.  "Multiple-Criteria  Decision  Making:  A 

Retrospective  Analysis."  IEEE  Transactions  on  Systems. 
Man,  and  Cybernetics.  SMC-15:  314  (May/June  1985). 

Hall,  A.D.  ”A  Three-Dimensional  Morphology  of  Systems 

Engineering.”  IEEE  Transactions  in  Systems,  Science,  and 
Cybernetics.  5:  156-lbO  (April  1969). 

A  Methodology  for  Systems  Engineering.  Princeton. 

MJ :  Van  Nostrand.  1962. 


A 


BIB-3 


Hard.  Colonel  Donald  G  .  .  Commander ,  L  S  AF  Satellite  Control 
Facility.  Personal  interview.  Air  lore e  Institute  ot 
Technology.  WPAFB  OH.  5  September  1985. 

Hazard.  Hap.  Li  .  S  ■  Space  Launch  Systems.  Report  No .  NSSAR- 
20-72-2,  Naval  Space  Systems  Activity.  Los  Angeies  CA  . 
'larch  1973  (AD  917  514). 

Heald.  Daniel  A.  ’’Economic  Benefits  of  the  OTV  Program.  ” 

Acta  Astronaut ica .  8:  1237-1249  (November-December 

1981  )  . 

Hill.  L.S.  ’Systems  Engineering  in  Perspective."  IEEE 

Transactions  on  Engineering  Management.  EM- 17:  124-131 

(November  1970). 

Hill.  Ph'lip  G  and  Carl  R.  Peterson.  Mechanics  an d 

Thet  in-Mi'  n  hi  !  cs  of  Propulsion.  Reading.  HA  Addison- 
Wesi>j'  t'uniishing  Company.  1970. 

Janson .  Lt  Col  Gotthard.  Director.  Development  and 

Operations  Manned  Spaceflight  Support  Group,  USAF  Det  2. 
Personal  interview.  NASA  Johnson  Space  Center,  Houston 
TX .  23  August  1985. 

Jensen.  Jorgen  et  at.  Design  Guide  to  Orbital  Flight.  New 
York:  McGraw-Hill  Book  Company,  Inc..  I9b2. 

Kamenetzky,  Ricardo  D.  ’’The  Relationship  between  the 
Analytic  Hierarchy  Process  and  the  Additive  Value 
Function.”  Decision  Sciences.  13:  702  (October  1982). 

Keeney.  Ralph  L.  and  Howard  Raiffa.  Decisions  with  Multiple 
Objectives-  Preferences  and  Value  Tradeoffs.  New  York: 
John  Wiley  and  Sons,  Inc..  1976. 

Kobayashi.  Hisashi.  Modeling  and  Analysis:  An  Introduction 
to  System  Performance  Evaluation  Methodology. 
Philippines-  Add i son-Wes 1 ey .  1978. 

Koelle.  Dietric  E.  ’’The  TRANSCOST  -  Model  for  Estimation  of 
Launch  Vehicle  Development.  Fabrication  and  Operations 
Cost,”  Advances  in  the  As t ronaut i ca 1  Sciences.  Volume 
49.  399-418.  linivelt.  Inc..  San  Diego,  1982. 

- .  ’’The  Impact  of  Launch  Vehicle  Type  and  Size  on 

Development  Cost."  Acta  As t ronaut i ca ,  8:  1195-1205 

(November-December  1981). 

Lemon.  George.  D  i  rector  .  D  i  rectora te  of  Logistics  Policy  and 
Programs.  Personal  interview.  HQ  SD.  Los  Angeles  AFS  CA . 
2  July  1985. 

Markowitz,  H.  M.  Portfolio  Selection.  New  York:  John  Wiley 
and  Sons.  Inc..  1959. 


BIB  1 


Martin.  Francis  F.  Computer  Modeling,  and  Simulation.  New 
York:  John  Wiley  and  Sons.  1 nc . .  19b8. 

Mylander.  W.  C.  et  a  1  .  A  Guide  tolSIMT-Version  4  The 
Computer  Program  Implementing  the  Sequential 
Unconstrained  Minimization  Technique  tor  Nonlinear 
Programm i ng ■  Research  Analysis  Corporation.  McLean  VA 
1971  (AD  731  391  )  . 

NASA.  1985  Long-Range  Program  plan.  NASA  Headquarters. 
Washington  D.C.:  October  1984a. 

NASA  Systems  Engineering  and  Integration  Space  Station 

Program  Office.  Space  Station  Reference  Configuration 
Description.  NASA  Lyndon  B.  Johnson  Space  Center. 
Houston  Texas:  August  1984b. 

NASA  Office  of  Aeronautics  and  Space  Technology.  Space 

Systems  Technology  Model:  Space  Technology  Trends  and 
Forcasts .  Vol  2  part  B,  NASA  Headquarters.  Washington 
D.C.  January  1984c. 

- •  Space  Systems  Technology  Model:  Analysis  of 

Technology  Needs  and  Future  Considerations.  Vol  3 .  N  A 
Headquarters,  Washington  D.C.  January  1984d. 

NASA  Office  of  Aeronautics  and  Space  Technology.  Space 
Systems  Technology  Model:  Executive  Summary.  NASA 
Headquarters.  Washington  D.C.  September  1981a. 

- .  Space  Systems  Technology  Model:  Mission  Systems 

Descriptions  and  Technology  Needs.  Vol  1  part  A.  NAS^ 
Headquarters.  Washington  D.C.  September  1981b. 

- •  Space  Systems  Technology  Model:  Instrument  Systems 

Description  and  Technology  Needs.  Vol  1  part  B.  NASA 
Headquarters.  Washington  D.C.  September  1981c. 

- •  Space  Systems  Technology  Model:  Opportunity 

Systems/Programs  and  Technology  .  Vol  3  .  NASA 
Headquarters,  Washington  D.C.  September  1981d. 

Pritsker,  Alan  A.  Introduction  to  Simulation  and  Slam  II 
(Second  Edition).  West  Lafeyette  Indiana:  Systems 
Publishing  Corporation.  1984. 

Rand  Corporation.  The.  Military  Equipment  Cost  Analysis. 
Santa  Monica.  CA :  The  Rand  Corporation.  June  1971 
(AD901477)  . 

Rehder  .  John  J.  and  Kathryn  E.  Wurster.  "Electric  vs 

Chemical  Propulsion  for  a  Large-Cargo  Orbit  Transfer 
Vehicle.  ”  J. Spacecraft .  1 b :  129-134  (May-June  1  d  7  9  ) 


II  I  II  5 


Saaty,  Thomas  L.  The  Analytical  Hierchv  Process.  Mew  York: 
McGraw-Hill  Book  Company.  Inc..  1980. 

Sage.  Andrew  P.  Methodology  for  Large  Scale  Systems.  New 
York:  McGraw-Hill  Book  Company.  1977. 

Space  Division.  Air  Force  Systems  Command.  Statement  of  Work 
for  the  Joint  DOD/NASA  Space  Transportation  Architecture 
Study .  Request  for  Proposal  No.  F0470 1 -85-R-0024 .  Los 
Angeles  CA .  31  May  1985. 

Stimpson.  W.  A.  Multiple  Attribute  Decision  Analysis  Model. 
ADA  111104  (Vol  1).  ADA  111105  (Vol  11  -  Programmer’s 
Guide).  DT1C.  Washington  D.C..  1983. 

^ndberg .  Lt  Col  Eric.  Personal  interview.  HU  SD.  Los 
Angeles  AFS  CA ,  3  July  1985. 

Tribus.  Myron.  Rational  Descriptions.  Decisions,  and 
Designs .  New  York:  Pergamon  Press.  Inc..  1969. 

Watkins.  Glen.  User  servicing  manager.  Space  Station. 

Telephone  interviews.  NAS"A  Headquarters.  Washington 
D.C..  9  August  through  12  August  1985. 

Wimberly.  Major  G.  V..  Chief , Advanced  Space  Communications 
Systems  Engineering.  Personal  interview.  HQ  SD .  Los 
Angeles  AFS  CA .  3  July  1985. 

Wismer,  David  A.  and  R.  Chattergy.  Introduction  to 

Nonlinear  Optimization:  A  Problem  Solving  Approach.  New 
York:  Elsevier  North  Holland.  Inc..  1978. 

Wittress.  Col  William.  Director  of  Plans.  Deputy  Commander 
for  Launch  and  Control  Systems.  Personal  interview.  HQ 
SD.  Los  Angeles  AFS  CA .  2  July  1985. 

Zeleny.  M.  ’’The  Theory  of  the  Displaced  Ideal.”  Multiple 
Criteria  Decision  Making:  Kyoto  1975.  (1975). 

Zersen.  Col  William  F.  H..  Assistant  for  Advanced  Launch 
Systems.  Deputy  Commander  for  Launch  and  Control 
Systems.  Personal  interview.  HQ  SD.  Los  Angeles  AFS  CA. 

2  July  1985. 


BIB-6 


Captain  Jeffrey  W.  Anderson  was  brn  on  1  July  1958  in 
Melbourne,  Florida.  He  graduated  from  high  school  in 
Beavercreek,  Ohio,  in  1976  and  attended  the  University  of 
Michigan  from  which  he  received  the  degree  of  Bachelor  of 
Science  in  Materials  and  Metallurgical  Engineering  in 
December  1980.  Upon  graduation,  he  received  a  commission  in 
the  USAF  through  the  ROTC  program.  He  was  called  to  active 
duty  in  February  1981.  He  served  as  integration,  test,  and 
evaluation  project  engineer  for  the  Sidewinder  Program 
Office  at  Naval  Air  Systems  Command  in  Washington,  D.  C. 
until  entering  the  School  o£  Engineering,  Air  Force  Insti¬ 
tute  of  Technology,  in  May  1984. 


Permanent  address:  2585  Brown  Bark  Dr. 


Fairborn,  Ohio  45324 


VITA 


Captian  Mark  S.  Gibson  was  born  on  8  Febuary  1954  in 
Milwaukee,  Wisconsin.  He  graduated  from  high  school  in 
Milwaukee,  Wisconsin,  in  1972.  In  1976  he  enlisted  in  the 
Air  Force.  In  1981  he  graduated  from  the  University  of 
Colorado,  with  a  Bachelor  o  Science  degree  in  Mechanical 
Engineering.  He  received  a  commission  through  the  Airman 
Educational  Commissioning  Program  and  Officer  Training 
School.  He  was  employed  as  an  astronaut ica 1  engineer  at 
Space  Division,  AFSC,  until  entering  the  School  of  Engineei — 
ing.  Air  Force  Institute  of  Technology,  in  May  1984. 


Permanent  Address:  2845  Loyola 

Riverside.  California 
92503 


* 

t 


VITA 


Captain  Delbert  B.  Langerock  was  born  on  16  November 
1954  in  Ei  Paso.  Texas.  He  graduated  from  high  school  in 
Marion,  South  Dakota  in  1973  and  entered  the  Air  Force  as  a 
Telecommunications  Control  Specialist.  Following  tours  of 
duty  at  RAF  Uxbridge  (Hillingdon),  England  and  NORAD  Combat 
Operations  Center,  Colorado,  he  was  selected  into  the  Airmen 
Education  and  Commissioning  Program  in  1978.  He  attended 
Iowa  State  University  from  which  he  received  the  degree  of 
Bachelor  of  Science  in  Industrial  Engineering  in  May  1981. 
Upon  graduation,  he  attended  USAF  Officer  Training  School 
and  received  a  commission  in  the  USAF  in  August  1981 .  He 
then  served  as  a  Minuteman  II  System  Technical  Engineer  in 
the  44th  Stragic  Missile  Wing,  Ellsworth  AFB,  South  Dakota, 
until  entering  the  School  of  Engineering,  Air  Force  Insti¬ 
tute  of  Technology,  in  May  1984. 

Permanent  address:  P.0.  Box  29 

Monroe,  South  Dakota 
57047 


VITA 


Captain  Richard  A.  Lieber  was  born  on  2b  September  1954 
in  Ponca  City,  Oklahoma.  He  graduated  from  high  school  in 
New  Orleans,  Louisiana,  in  1972  and  attended  Oklahoma  State 
University  from  which  he  received  the  deree  of  Bachelor  of 
Science  in  Electrical  Engineering  in  December  197b.  Upon 
graduation,  he  received  a  commission  in  the  USAF  through  the 
ROTC  program.  He  was  employed  as  a  staff  manager  and  pro¬ 
duction  supervisor  for  the  Procter  &  Gamble  Manufacturing 
Company,  St.  Louis.  Missouri,  until  called  to  active  duty  in 
January  1978.  He  completed  pilot  training  and  received  his 
wings  in  December  1978.  He  then  served  as  a  copilot  and 
aircraft  commander  in  the  325  Bombardment  Squadron.  Fair- 
child  AFB.  Washington,  until  entering  the  School  of 
Engineering,  Air  Force  Institute  of  Technology,  in  May  1984. 


Permanent  address:  7bl8  Teal  Run  Drive 


Houston.  Texas  77071 


VITA 


Captain  Michael  A.  Palmer  was  born  on  30  August  1959  in 
Highland  Park,  Michigan.  He  graduated  from  high  school  in 
Lathrup  Village,  Michigan,  in  1977  and  attended  Oakland 
University  from  which  he  received  the  degree  of  Bachelor  of 
Science  in  Mechanical  Engineering  in  April,  1981.  Following 
graduation,  he  received  a  commission  in  the  USAF  after 
attending  OTS .  He  then  served  as  Mission  Planning  Project 
Engineer  for  the  6595  Shuttle  Test  Group,  Vandenberg  AFB, 
California,  in  preparation  for  the  first  Space  Shuttle 
launch  from  the  Western  Test  Range.  He  then  entered  the 
School  of  Engineering,  Air  Force  Institute  of  Technology  in 
May,  1984. 


Permanenet  Address:  2218  Kipling 

Sterling  Heights,  Michigan 


48078 


Captain  Michael  W.  Peltzer  was  born  23  February  1959  in 
North  Hollywood,  California.  He  graduated  from  high  school 
in  Panorama  City,  California,  in  1976  and  attended  the 
University  of  California,  at  Los  Angeles  (UCLA)  from  which 
he  received  the  degree  of  Bachelor  of  Science  in  Mechanical 
Engineering  in  December  1980.  Upon  graduation,  he  received 
a  commission  in  the  USAF  through  the  ROTC  program.  He 
entered  active  duty  in  January  1981  and  was  assigned  to  the 
Directorate  of  Materiel  Management,  Item  Management  Division 
of  the  Ogden  Air  Logistics  Center,  Hill  AFB,  Utah,  where  he 
served  as  a  landing  gear  systems  engineer.  He  entered  the 
School  of  Engineering,  Air  Force  Institute  of  Technology,  in 
May  1984. 

Permanent  Address:  8160  Wakefield  Avenue 

Panorama  City,  California 


SECURITY  CLASS  *  'CATION  Qf  This  page 


REPORT  DOCUMENTATION  PAGE 


report  stcu^iT v  classification 

UNCLASSIFIED 


lb  RESTRICTIVE  MARKINGS 


2«  security  Classification  authority 


2b  OEClASSiFiCAT.ON-OOWNGRAOiNG  SCHEDULE 


3  OISTRIBUTION/AVAILABILITY  of  report 

Approved  for  public  release;  distribution 
unlimited. 


*  PERFORMING  ORGANIZATION  REPORT  NUMBS  R(S) 

AFIT/GSE/ENY/85D 


5.  MONITORING  ORGANIZATION  REPORT  NUM8ERIS) 


6a  NAME  OF  PERFORMING  ORGANIZATION 

School  of  Engineering 
Air  Force  Institute  of  Technology  AFIT/ENY 


|SB  OFFICE  SYMBOL 
(If  applicable ) 


7a.  NAME  OF  MONITORING  ORGANIZATION 


6c  ADDRESS  (City  State  and  ZIP  Code > 

Wright-Patterson  AFB  Ohio,  45433 


7b.  ADORESS  (City.  State  and  ZIP  Code ) 


la  NAME  OF  FUNOING/SPONSORING 

ORGANIZATION 

Office  of  manned  spaceflight 
Office  of  plans _ 


1 8b.  OFFICE  SYMBOL 

(if  applicable > 

SPAM,  SD/XR 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


Be  ADORESS  'City.  State  and  ZIP  Code) 


10.  SOURCE  OF  FUNDING  NOS. 


P.O.  Box  92960,  Worldway  Postal  Center 

PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

Los  Angeles,  CA  90009-2960 

element  no 

NO. 

NO. 

NO 

11  TITLE  'Include  Security  C loan ficat tom 

See  Box  19 

35171F 

662861 

088000 

00000 

12  PERSONAL  AUTHORISI 

"ee  Box  19 


JO.  TYPE  OF  REPORT 

13b  time  COVERED 

1«  DATE  OF  REPORT  lYr  .  Mo..  Day/ 

15.  PAGE  COUNT  | 

MS  Thesis 

FROM  TO 

1985  December 

111  | 

16  supplementary  notation 

Thesis  is  in  three  volumes:  Volume  I  -  Executive  Surmary,  Volume  II  -  Final  Report,  and 

Volume  III  -  Appendices. 


1 1 7  COSATI  CODES  | 

FIELD 

GROUP 

SUB  GR 

22 

01 

18.  SUBJECT  TERMS  (Continue  on  reverie  if  neceuary  and  identify  by  block  number) 

SATELLITES,  REPAIR,  SPACE  MAINTENANCE 


19  ABSTRACT  iContmue  on  reverie  if  neceiaary  and  identify  by  block  number i 

Title:  A  METHODOLOGY  FOR  SELECTION  OF  A  SATELLITE  SERVICING  ARCHITECTURE 
VOLUME  III  -  APPENDICES 

Personal  Authors:  Jeffrey  W.  Anderson,  Capt  USAF 
Mark  S.  Gibson,  Capt  USAF 
Delbert  B.  Langerock,  Capt  USAF 
Richard  A.  Lieber,  Capt  USAF 
Michael  A.  Palmer,  Capt  USAF 
Michael  W.  Peltzer,  Capt  USAF 


MS*  u«  wb  m-j f. 

e .  woiv^  f/» 

lor  R«i*ouh  and  Profowlona!  *D«T«!epaaa| 
Air  Forco  Instil  M  Twch:i  <y 
Wrtflht-Potterfc.  n  <  •  U  U/i 


Thesis  Chairman:  Stuart  C.  Kramer,  Capt  USAF 

Assistant  Professor  of  Aerospace  Engineering 


OISTRIBUTION/AVAILABILITY  of  abstract 

Iunclassifieo/unlimiteoXS  same  as  rpt  □  otic  users  □ 


21  abstract  security  classification 

UNCLASSIFIED 


I  23*.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Stuart  C.  Kramer,  Capt  USAF 


22b  TELEPHONE  NUMBER 


22c  OFFICE  SYMBOL 

AFIT/ENY 


UNCLASSIFIED 


IfCUHiTV  CLASSIFICATION  OF  THIS  FAGS 


VOLUME  III  -  APPENDICES 

A  two-phase  methodology  for  selecting  an  optimal  military 
satellite  servicing  system  is  developed  using  the  systems  engin¬ 
eering  approach.  This  methodology  is  used  to  evaluate  several 
alternative  systems  at  varying  levels  of  detail.  The  candidate 
systems  cure  ocitposed  of  lew-G  launchers,  high-G  launchers,  orbital 
servicing  vehicles,  and  space  bases.  An  optinal  realization  is 
then  derived  for  a  system  of  lcwH3  launchers  and  orbited,  servic¬ 
ing  vehicles.  In  the  first  phase  of  the  approach,  vector  optimiz¬ 
ation  techniques  are  used  to  vary  the  states  of  a  model  to  obtain 
a  set  of  optimal  solutions.  The  second  phase  embodies  the  decision 
maker's  preferences  in  a  value  system  to  enable  preference  ranking 
of  the  optinal  solutions  in  the  non-daminatad  solution  set.  This 
methodology,  as  presented,  can  be  applied  to  any  ccrrplex  problem 
with  multiple  conflicting  objectives.  It  is  designed  for  use  by 
an  engineering  organization  supporting  a  senior- level  decision 
naker. 

The  report  is  in  three  volutes.  The  Executive  Suntary  (Volute  I) 
is  a  cursory  review  of  the  study  and  is  meant  to  be  self-contained. 
The  Final  Report  (Volume  II)  and  the  Appendices  (Volute  III)  are  more 
detailed  and  should  be  read  together  for  ocnpleteness . 


TOOJSSSTFIHT 


SICUNITV  CLASSIFICATION  OF  THIS  FAOF 


DATE 
ILME 


